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Description 

[0001] The present invention relates to a spinning strip aperture radiometer sensor system comprising a telescope 
comprised of a rotating strip aperture that rotates around an optical axis and that produces temporarily sequential 
5 images of an image scene, a two-dimensional detector array for detecting images formed at the focal plane of the 
telescope, and a signal processor coupled to the detector array for recording a plurality of image frames of the image 
scene as the strip aperture rotates around the optical axis of the telescope, and for synthesizing images from the 
plurality of image frames. 

[0002] The present invention relates further to an object detection method for use in an imaging system comprising 
10 a spinning strip aperture telescope that rotates around an optical axis, a two-dimensional detector array for detecting 

images located at a focal plane of the telescope, and a signal processor for processing the detected images, wherein 

the object detection method comprises the following steps: 

recording a plurality of image frames of an image scene imaged by the telescope as the strip aperture rotates 

around the optical axis of the telescope, and synthesizing images from the plurality of recorded image frames. 
15 [0003] A spinning strip aperture radiometer sensor system and such an object detection method are known from 

US-A-5,243,351. 

[0004] The present invention relates generally to spinning aperture radiometers and methods, and more particularly 
to spinning strip (partial) aperture imaging radiometers that synthesize equivalent full circular aperture images and 
detect moving objects from a plurality of rotating strip aperture image, or synthesized image measurements. 

20 [0005] To provide high resolution images from space-based platforms, for example, conventional sensor architec- 
tures use active control of large, heavy, deployable optical systems. Depending upon the mission requirements and 
size of the primary mirror, the active control can range from periodic piston and tilt control of primary mirror segments 
to piston, tilt, figure, and alignment control of all optical elements of the sensor. Full aperture systems with the same 
resolution as the present invention have a great deal of light gathering capability because of the relatively large aperture 

25 areas. However, to place multi-meter diameter apertures into orbit, full aperture systems competing with the present 
invention require: segmented optical surfaces and folded support structures, if the optical system diameters are larger 
than the launch vehicle's fairing; complex and potentially high bandwidth adaptive optical techniques, if thin deformable 
mirrors are used to save weight; and complex piston and pupil matching control, if implemented as a phased array. 
Therefore, full aperture systems are relatively heavy and have high technical risk when compared to the present in- 

30 vention. 

[0006] The processors and methods described herein relate to and provide for an improvement to the system capa- 
bilities of the closest prior art. referred to herein as SpinAp. The SpinAp system and processing method is described 
in the above-mentioned U.S. Patent No. 5,243,351, entitled "Full Aperture Image synthesis Using Rotating Strip Ap- 
erture Image Measurements", which is assigned to the assignee of the present invention. 

35 [0007] SpinAp is a sensor system and data processing method that is capable of synthesizing images having reso- 
lution equivalent to a full circular aperture. The equivalent SpinAp synthesized full circular aperture has a diameter 
equal to the largest correlation length associated with the strip aperture's geometry. To accurately synthesize an equiv- 
alent full aperture image, SpinAp acquires imaging data using the rotating sensor's focal plane detectors at time intervals 
appropriate for completely measuring the optical system passband of the equivalent full circular aperture. The images 

^0 are processed by methods described in U.S. Patent No. 5,243,351. 

[0008] The commonality of the approaches between the SpinAp system and the system of the present invention 
(referred to as SpinAp moving object detection) arises from the use of temporally registered strip aperture measure- 
ments to synthesize images having resolution equivalent to a full circular aperture. The equivalent SpinAp synthesized 
full circular aperture has a diameter equal to the largest correlation length associated with the strip aperture's geometry. 

45 [0009] To accurately synthesize an equivalent full aperture image, SpinAp acquires imaging data using the rotating 
sensor's focal plane detectors at time intervals appropriate for completely measuring the optical system passband of 
the equivalent full circular aperture. The images are processed by methods described in U.S. Patent No. 5,243,351 . 
[0010] However, a moving object in the instantaneous field of view of any of the individual frames will result in a 
motion smear signature in the total synthesized image. Detection of the moving object is required to establish object 

50 phenomenology, such as, position, size, trajectory, velocity, acceleration, and point of origin. In addition, the detection 
process is a precursor to applying processors and methods to the acquired data that permit stationary image synthesis 
of the moving objects. Therefore, a hierarchy of moving object detection processors and methods is described herein. 
The hierarchy includes spatial, temporal, spatial frequency, and temporal frequency domain detection processors, 
which may also incorporate multi-spectral background rejection techniques. The present system makes use of the 

55 available measurements to detect moving objects in the field of view of the synthesized images, in the field of view of 
any pair of SpinAp individual frames, in any one of the SpinAp individual frames, or any combination. The same methods 
may also be applied to transient event detection. 

[0011] Therefore, it is an objective of the present invention, to provide an improved SpinAp system and method. 
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[0012] This objective is achieved by the radiometer sensor system, mentioned at the outset, wherein the signal 
processor is programmed for differencing pairs of images or frames, and wherein a matched filter is provided, so as 
to detect moving objects found in the synthesized images. 

[001 3] The above objective is further achieved by the object detection method mentioned at the outset, wherein the 
5 method comprises the steps of differencing pairs of imags or frames, and processing the resulting data through a 
matched fitter, so as to detect moving objects found in the synthesized images. 

[0014] It is an advantage to provide a modified SpinAp image processing method that results in lower payload weight 
for a given effective synthesized aperture size, while providing capability to detect moving objects in the sensor's field 
of view. Furthermore, it is advantageous to have a system having image processing and moving object detection 
10 methods that provide a lower risk, lower cost, lighter weight, and simpler fabrication deployment alternatives to de- 
ploying complex, large full circular apertures (or phased array telescopes), requiring complex adaptive optical systems 
for space based imaging applications. 

[0015] In general, to meet the above and other objectives, one embodiment of the present invention provides for a 
spinning strip radiometer imaging and detection system that is capable of detecting moving objects. The system in- 

*5 eludes a rotating strip aperture telescope that produces temporally contiguous or sequential images. The rotating strip 
aperture telescope typically comprises a rotating strip aperture primary reflector and a secondary reflector. A two di- 
mensional detector array is provided to detect images located in the focal plane of the telescope. A rotation compen- 
sation device is employed to prevent rotational smear during the integration time of the detectors of the array. A signal 
processor is provided for recording a plurality of image frames of a target scene imaged by the telescope as the strip 

20 aperture rotates around the optical axis of the telescope, for synthesizing the equivalent full circular aperture image, 
and for detecting moving objects in'the individual strip aperture frames, any pair of frames, the synthesized image, or 
any combination of individual frames and synthesized images. 

[0016] The present invention thus provides for a spinning strip (partial) aperture imaging system that synthesizes 
equivalent full circular aperture images, and detects moving objects in the sensors field of view from a plurality of 
25 rotating strip aperture image measurements, while compensating for random, and/or systematic line of sight errors 
between individual strip aperture images. The present invention thus provides improved total system performance 
when compared to the SpinAp system of U.S. Patent No. 5,243.351 by providing means to synthesize an equivalent 
full circular image, and simultaneously providing the capability to detect moving objects. 

[0017] Example embodiments of SpinAp moving object detection methods are described below. Since each embod- 
30 iment of the radiometer system, image synthesis, and moving object detection processing method depends upon spe- 
cific mission requirements and engineering tradeoffs, the radiometer system; image synthesis, and moving object 
detection method incorporates means to compensate for random, and/or systematic line of sight drift between frames, 
and a priori and a posteriori known error sources, such as. non-isoplanatic optical system point spread functions, field 
point independent image smear due to image motion and finite electronic bandwidth of the detectors, and field point 
35 dependent image smear caused by uncompensated rotational motion of the image. 

[0018] Embodiments of the image synthesis and detection process performed by the present invention summarized 
above are as follows. As the spinning strip aperture rotates around the telescope's optical axis the following occurs. 
The rotation compensation device counter-rotates during the detector integration time, thereby providing a temporally 
stationary image. An image frame is recorded and saved. If a rotating (relative to the scene) detector array architecture 
40 has been selected, the acquired frame is coordinate transformed and interpolated to a reference grid of the synthesized 
image. The data is Fourier transformed and stored. A new frame is recorded and saved. 

[0019] An estimate of the frame-to-frame misregistration of the recorded data due to random line of sight errors is 
obtained. The strip aperture images, or the Fourier transforms, are corrected for their line of sight errors and are stored. 
The preceding steps are sequentially repeated for each strip aperture image frame, or the frames are sequentially 

45 acquired and stored, and then global estimates of the line of sight errors are determined and used to register the 
frames. Once any pair, or combination of pairs, of the individual strip aperture images have been corrected for the line 
of sight errors, the SpinAp moving object detection process may begin. In addition, once a sufficient number (as de- 
scribed in U.S. Patent No. 5,243,351 ) of strip aperture image measurements have been acquired the image synthesis 
process may begin. The SpinAp moving object detection method may be performed on two, or more, synthesized 

so images. 

[0020] An embodiment of the SpinAp moving object detection method using two or more synthesized images acquired 
by observing the same scene spatially registers the synthesized images (if required due to line of sight errors between 
measurement sets), differences all possible distinct pairs of synthesized images, sums the differenced images, and 
passes the resulting differenced and summed data through a matched filter array and then a comparator. A moving 
55 object detection is declared when any of the outputs from the matched filter array exceed the predetermined thresholds 
of the comparator. 

[0021] Another embodiment of the SpinAp moving object detection method using two or more individual line of sight 
corrected frames resulting from the observation of the same scene applies a system transfer function (STF) equalization 
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filter to the individual frames to be differenced, differences any or all distinct pairs of equalized frames, and applies a 
non-common passband rejection filter to the resulting differenced data. The resulting filtered and differenced data are 
combined by addition to form the total summed spectrum of filtered differences. Trie summed spectra are inverse 
transformed, thereby providing a moving object signature in the spatial domain, which is transferred to a noise optimal 
5 matched filter array and a comparator. A moving object detection is declared when any of the outputs from the matched 
filter array exceed the predetermined thresholds of the comparator. 

[0022] The noise optimal matched filter detection processor incorporates properties corresponding to the physical 
phenomenology of potential objects, such as, shape, size, velocity, acceleration, background to object contrast ratio, 
signal to noise ratio, number of objects, number of track crossings, and the like. The matched filter detection processor 

10 also incorporates a priori knowledge of the measurement noise statistics propagated through the equalization, differ- 
encing, and filtering process, as well as the unique signal measurement characteristics of the SpinAp sensor system. 
[0023] The equalization, differencing, and filtering operations associated with the present invention eliminate all tem- 
porally stationary components in a common overlap passband of the individual frame images, as well as the synthesized 
images, leaving only a noisy signature of smeared moving objects. The performance of the moving object detection 

15 techniques is scenario dependent and is degraded by individual frame measurement noise propagated through the 
equalization, differencing, and filtering operations, and background leakage due to residual frame to frame, and/or 
synthesized image to synthesized image registration error. 

[0024] Alternative embodiments of spatial (and spatial frequency) methods that are described herein correspond to 
applying optimal matched filters to any of the differenced frames, and combining the outputs for all possible frame 
20 differences. 

[0025] While the previous embodiments are typical of SpinAp temporal domain detection methods, temporal fre- 
quency domain processes may also be employed. An embodiment of a temporal frequency domain processor utilizes 
a temporal sequence of the spatially registered individual SpinAp strip aperture frames obtained by observations of 
the same scene. The spatially registered frames are spatially Fourier transformed to a common spatial frequency grid. 

25 The common spatial frequency grid is associated with the equivalent full circular aperture SpinAp synthesized spatial 
frequency passband. Any or all possible spatially registered and Fourier transformed individual frame images are tem- 
porally Fourier transformed to produce a data set associating all temporal frequencies with each spatial frequency. To 
obtain the best possible detection performance, the temporal and spatial frequency transformed data set is transferred 
to a noise optimal matched filter detection processor. The matched filter detection processor incorporates a priori 

30 knowledge of the object phenomenology, signal and noise measurement statistics, as well as the unique measurement 
characteristics of the SpinAp sensor system. A moving object detection is declared when any of the outputs from the 
matched filter array exceed predetermined thresholds. 

[0026] By eliminating the zero temporal frequency, the matched fitter implicitly eliminates all temporally stationary 
components in the transformed time history of the individual frames, thereby leaving only the temporal frequency and 
35 spatial frequency signature of the moving objects. The performance of the moving object detection technique is scenario 
dependent and is degraded by individual frame measurement noise propagation, and background leakage due to 
residual frame to frame registration error. 

[0027] If observational parameters permit, acceptable performance may be achieved for less than optimal approach- 
es, whereby the optimal matched filter array and comparator can be replaced by template matched filters and a com- 
40 parator. streak detectors and a comparator, or a simple single pixel or multiple pixel threshold exceedence detection 
processor, In this context the word pixel refers to the spatial, or the spatial frequency element of the two dimensional 
discrete representation of the image, or its transform. 

[0028] A hierarchy of moving object detection processors and methods is provided by the present invention. The 
hierarchy includes spatial, temporal, spatial frequency, and temporal frequency domain detection processors, as well 
45 as, multi-spectral background rejection techniques. Selection of the appropriate processing procedure and method 
depends upon the scenario, and the effective signal to noise ratio characteristics of the moving object. The same 
methods may be applied to transient event detection. 

[0029] * Given knowledge of the strip aperture's spatial response function, the spatial response function of the detec- 
tors of the array, the noise statistics, and the temporal registrations of each of the recorded strip aperture images 
so permits moving object detection by the sensor system and processing methods. 

[0030] The various features and advantages of the present invention may be more readily understood with reference 
to the following detailed description taken in conjunction with the accompanying drawings, wherein like reference nu- 
merals designate like structural elements, and in which: 

55 Fig. 1 illustrates a spinning aperture imaging radiometer system employing improved data processing methods in 

accordance with the principles of the present invention; 

Fig. 2 illustrates a top level embodiment of a SpinAp moving object detection processing method in accordance 
with the principles of the present invention utilizing a plurality of spatially registered SpinAp synthesized images; 
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Fig. 3 illustrates a top level embodiment of a SpinAp moving object detection processing method in accordance 
with the principles of the present invention utilizing a plurality of spatially registered, system transfer function equal- 
ized individual SpinAp transformed frames; 

Fig. 4 illustrates a top level embodiment of a SpinAp moving object detection processing method in accordance 
5 with the principles of the present invention utilizing a plurality of spatially registered, temporal and spatial frequency 

transforms of SpinAp individual frame images; and 

Fig. 5 illustrates a top level embodiment of a SpinAp moving object detection processing method in accordance 
with the principles of the present invention utilizing a plurality of spatially registered, system transfer function equal- 
ized, individual SpinAp transformed frames, and sequential difference frame matched filtering. 

10 

[0031] Referring to the drawing figures, Fig. 1 illustrates a spinning aperture imaging radiometer system 10 and data 
processing methods 20 (corresponding to the system described in U.S. Patent No. 5,243,351) that is appropriately 
modified using the principles of the present invention. The spinning aperture imaging radiometer system 10 and data 
processing method 20 provides moving object detection information, and provide the capability to generate equivalent 

15 full circular imagery, while removing line of sight jitter and frame registration errors. Therefore, the system 10 and data 
processing methods 20 described herein provide enhancements to the SpinAp system of U.S. Patent No. 5,243,351 . 
[0032] The spinning aperture imaging radiometer system 10 comprises a rotating strip aperture telescope 11 that 
comprises a primary 12a and a secondary reflector 12b. A tertiary reflector (not shown) may be employed in the tele- 
scope 11 under certain circumstances. The system 1 0 is shown in the form of a satellite comprising a stationary section 

20 13 having earth pointing antenna 13a. The telescope 11 is disposed on a platform 14, to which the stationary section 
1 3 is also coupled. The spinning aperture imaging radiometer system 1 0 is adapted to record a number of image frames 
of a target scene imaged by the telescope 11 as the primary mirror 1 2a (comprising a strip aperture 12) rotates around 
the telescope's optical axis. 

[0033] A line of sight stabilization mirror 1 5 and an image derotation device 16 are disposed along the optical path 
25 of the telescope 1 1 that are adapted to stabilize and derotate the image prior to its sensing by a detector array 1 7. The 
derotation device 16 counter rotates the image during the integration time of detectors comprising the detector array 
17, under control of a rotation compensation controller 19, thereby providing a stationary image. The line of sight 
stabilization mirror 15 is used by a line of sight control system (such as may be provided by a signal processor 18 or 
other dedicated control system) to remove high bandwidth line of sight errors, as well as line of sight errors due to 
30 orbital dynamics of the system 10. 

[0034] The target scene is imaged onto the detector array 17 located at the focal plane of the telescope 11 that is 
coupled to the signal processor 18 that is adapted to process the image frames. Individual image frames are processed 
and combined in the signal processor 18 to synthesize equivalent full circular aperture images. In the present enhanced 
moving object detection SpinAp system 1 0, the signal processor 1 8 (or an independent processor) is modified according 
35 to the methods of the present invention to provide for moving object detection. Detailed descriptions of the modified 
processing methods are illustrated with reference to Figs. 2-5. 

[0035] Figs. 2-5 illustrate a variety of different embodiments of the sensor system 10 and processing methods 20, 
20-2, 20-3, 20-4 in accordance with the principles of the present invention employed by the spinning aperture radiometer 
system. A particular embodiment is selected based upon the observational scenario and timelines, the amount of 

*o available a priori knowledge, and the available computational throughput of the processing chain. 

[0036] Fig. 2 illustrates a top level diagram of a simple SpinAp moving object detection processing method 20 in 
accordance with the principles of the present invention employed using an appropriately modified spinning aperture 
radiometer system and processing methods described in U.S. Patent No. 5,243,351 . The detection method 20 shown 
in Fig. 2 employs two or more SpinAp strip aperture image sets that observe the same scene to generate equivalent 

« synthesized full circular aperture images 30 of the scene in accordance with U.S. Patent No. 5,243,351 . In the presence 
of line of sight registration errors, each of the synthesized images 30 are optimally registered spatially by methods 
described in U.S. Patent No. 5,243,351. 

[0037] Any, or all, possible distinct pairs of spatially registered synthesized images 30 are differenced 40 and sub- 
sequently added 35 and stored, thereby generating a total sum 45 of differenced synthesized images. If the number 
50 of distinct synthesized images 30 is N, then the total number of possible distinct differenced synthesized image pairs 
is given by N(N-1)/2. 

[0038] To obtain the best possible detection performance, the resulting sum 45 of differenced synthesized images 
is transferred to a noise optimal matched filter detection processor 50. The matched filter detection processor 50 in- 
corporates an array of matched filters 55 having properties corresponding to the physical phenomenology of potential 
55 objects, such as, shape, size, velocity, acceleration, background to object contrast ratio, signal to noise ratio number 
of objects, number of track crossings, and the like. The matched filter detection processor 50 also incorporates a priori 
knowledge of measurement noise statistics propagated through the synthesis and differencing process, as well as, 
knowledge concerning the unique signal measurement characteristics of the SpinAp sensor system 1 0. Each of the 
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matched filters 55 of the detection processor 50 is applied to every pixel location of the total sum 45 of differenced 
synthesized images. The output of each matched filter 55 is transferred to a comparator 60 with thresholds set to a 
predetermined level of each matched filter 55. A moving object detection is declared when any of the outputs from the 
matched filters 55 exceed the predetermined thresholds of the comparator 60. The threshold may be predetermined 
5 by a priori assumptions concerning the signal and noise statistics, or by a posteriori assumptions concerning the signal 
and noise statistics based upon the individual frame or set of SpinAp frames. 

[0039] The differencing operation eliminates all temporally stationary components of the synthesized image, leaving 
only the signature of smeared synthesized moving objects. The performance of the moving object detection processing 
method 20 is scenario dependent and is degraded by individual frame measurement noise propagated through the 
10 synthesis and differencing operations, and background leakage due to residual frame to frame registration error. If 
observation parameters permit, acceptable performance may be achieved for less than optimal approaches, whereby 
the optimal matched filters 55 and the comparator 60 may be replaced by template matched filters and the comparator 
60, streak detectors and the comparator 60, or a simple single pixel or multiple pixel threshold exceedence detection 
processor. 

15 [0040] Fig. 3 illustrates a top level embodiments a SpinAp moving object detection processing method 20-2 in 
accordance with the principles of the present invention employed using an appropriately modified spinning aperture 
radiometer system and processing methods described in U.S. Patent No. 5,243,351 . The detection method 20-2 shown 
in Fig. 3 employs two or more SpinAp strip aperture images 30-2 obtained by observing the same scene, and generated 
using methods described in U.S. Patent No. 5,243,351 . In the presence of line of sight registration errors, each of the 

20 individual SpinAp images 30-2 are optimally registered spatially as described in U.S. Patent No. 5,243,351 , and Fourier 
transformed to a common spatial frequency coordinate system and grid. Typically, the common spatial frequency grid 
is associated with an equivalent full circular aperture SpinAp synthesized spatial frequency passband. 
[0041] Any or all possible distinct pairs of spatially registered and Fourier transformed individual frame images 30-2 
are multiplied by spatial frequency filters 31-2 designed to compensate for the temporally varying system transfer 

25 functions associated with each individual frame and resulting from the rotation of the SpinAp sensor system 1 0. The 
system transfer function compensation may be implemented by one of several methods, including: equalizing each 
frame to a target full circular aperture or other reference system transfer function, equalizing one frame to the system 
transfer function of the frame to be differenced, or deconvolving the system transfer function from both images to be 
differenced. For ideal system transfer functions, and if and only if only frames acquired n radians apart are to be 

30 differenced, no system transfer function equalization is required. 

[0042] For a template signal detection approach, SpinAp moving object detection methods 20-2 based upon any 
one of the many possible equalization procedures have unique noise performance characterized by the particular 
method's receiver operating characteristics (ROC). The differences in ROC for signal template approaches are dom- 
inated by different and uncompensated levels of noise "boosting" or "de-boosting" generated by the equalization proc- 

35 ess. However, an appropriately designed noise optimal matched filter approach compensates for the noise character- 
istics in the equalization process. 

[0043] The equalized individual frames are differenced 40-2, and spatial frequency filtered 41-2 to eliminate non- 
common information outside an overlap passband associated with the two frames. The spatial frequency filtered frame 
difference is added 35-2 to previous filtered frame differences and stored, thereby generating a total sum 35-2 of 
<o differenced spatial frequency filtered Fourier transformed images. If the number of distinct individual frame images is 
N, then the total number of possible distinct differenced individual frame image pairs contributing to the sum may be 
as large as N(N-1)/2. 

[0044] To obtain the best possible detection performance, the resulting sum 35-2 of differenced spatial frequency 
filtered transformed images is inverse Fourier transformed 42-2 to obtain a spatial domain image 45-2, which is trans- 

45 ferred to a noise optimal matched filter detection processor 50-2. The matched filter detection processor 50-2 incor- 
porates an array of matched filters 55-2 having properties corresponding to the physical phenomenology of potential 
objects, such as, shape, size, velocity, acceleration, background to object contrast ratio, signal to noise ratio, number 
of objects, number of track crossings, and the like. The matched filter detection processor 50-2 also incorporates a 
priori knowledge of the measurement noise statistics propagated through the equalization, differencing, and filtering 

so process, as well as the unique signal measurement characteristics of the SpinAp sensor system 10. Each of the 
matched filters 55-2 of the detection processor 50-2 is applied to every pixel location of the total sum 35-2 of differenced 
filtered images. The output of each matched filter 55-2 is transferred to a comparator 60-2 with thresholds set to a 
predetermined level for each matched filter 55-2. A moving object detection is declared when any of the outputs from 
the array of matched filters 55-2 exceed the predetermined thresholds of the comparator 60-2. The threshold may be 

55 predetermined by a priori assumptions concerning the signal and noise statistics, or the a posteriori assumptions 
concerning the signal and noise statistics based upon the individual frame or set of SpinAp frames. 
[0045] The equalization, differencing, and filtering operations eliminate all temporally stationary components in the 
common overlap passband of the individual frame images, leaving only the signature of smeared moving objects. The 
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performance of the moving object detection processing method 20-2 is scenario dependent and is degraded by indi- 
vidual frame measurement noise propagated through the equalization, differencing, and filtering operations, and back- 
ground leakage due to residual frame to frame registration error. If observational parameters permit, acceptable per- 
formance may be achieved for less than optimal approaches, whereby the optimal array of matched filters 55-2 and 

5 comparator 60-2 may be replaced by template matched filters and the comparator 60-2, streak detectors and the 
comparator 60-2, or a simple single pixel or multiple pixel threshold exceedence detection processor. 
[0046] Fig. 4 illustrates a top level embodiment of a SpinAp moving object detection processing method 20-3 in 
accordance with the principles of the present invention employing an appropriately modified spinning aperture radiom- 
eter system and processing methods described in U.S. Patent No. 5,243,351 . The detection method 20-3 shown in 

10 Fig. 4 employs a temporal sequence of individual SpinAp strip aperture images 50-3 or sub-images obtained by ob- 
serving the same scene, generated using methods described in U.S. Patent No. 5,243,351 . In the presence of line of 
sight registration errors, each of the individual SpinAp images 50-3 are optimally registered spatially by techniques 
described in U.S. Patent No. 5,243,351, and Fourier transformed 51-3 to a common spatial frequency coordinate 
system and grid, thereby producing Fourier transformed images 52-3. Typically, the common spatial frequency grid is 

15 associated with an equivalent full circular aperture SpinAp synthesized spatial frequency passband. Any or all possible 
spatially registered and Fourier transformed individual frame images 52-3 are temporally Fourier transformed 53- 3 to 
produce a data set 54-3 associating all temporal frequencies with each spatial frequency. 

[0047] To obtain the best possible detection performance, the temporal frequency and spatial frequency transformed 
data set is transferred to a noise optimal matched filter detection processor 55-3. The matched filter detection processor 

20 55-3 incorporates an array of matched filters 56-3 having properties corresponding to the physical phenomenology of 
potential objects, such as, shape, size, velocity, acceleration, background to object contrast ratio, signal to noise ratio, 
number of objects, number of track crossings, and the like. The matched filter detection processor 55-3 also incorpo- 
rates a priori knowledge of the variation of the system transfer function and other unique measurement characteristics 
associated with the SpinAp sensor system 10, and the measurement signal and noise statistics propagated through 

25 the measurement and transforming processes. Each matched filter 56-3 of the detection processor 55-3 is applied to 
the set 54-3 of temporally and spatially transformed data. The output of each matched filter 56-3 is transferred to a 
comparator 60-3 with thresholds set to a predetermined level for each matched filter 56-3. A moving object detection 
is declared when any of the outputs from the array of matched filters 56-3 exceed the predetermined thresholds of the 
comparator 60-3. The threshold may be predetermined by a priori assumptions concerning the signal and noise sta- 

30 tistics, and/or by a posteriori assumptions concerning the signal and noise statistics based upon the individual frame 
or set of SpinAp frames. 

[0048] The temporal frequency matched filter 56-3 implicitly eliminates all temporally stationary components in the 
individual frame images, leaving only the temporal and spatial frequency signature of smeared moving objects. The 
performance of the moving object detection processing method 20-3 is scenario dependent and is degraded by indi- 
35 vidual frame measurement noise propagation, and background leakage due to residual frame to frame registration 
error. If observational parameters permit, acceptable performance may be achieved for less than optimal approaches, 
whereby the optimal array of matched filters 56-3 and comparator 60-3 may be replaced by template matched filters 
and the comparator 60-3, or a simple single frequency "pixel" or multiple frequency "pixel" threshold exceedence de- 
tection processor. 

40 [0049] Fig. 5 illustrates a top level embodiment of a SpinAp moving object detection processing method 20-4 in 
accordance with the principles of the present invention employing an appropriately modified spinning aperture radiom- 
eter system and processing methods described in U.S. Patent No. 5,243,351. The detection method 20-4 of Fig. 5 
employs two or more SpinAp strip aperture images 30-4 obtained by observing the same scene, generated in accord- 
ance with methods described in U.S. Patent No. 5,243,351 . In the presence of line of sight registration errors, each of 

45 the individual SpinAp images are optimally registered spatially by methods described in U.S. Patent No. 5,243,351, 
and Fourier transformed to a common spatial frequency coordinate system and grid. Typically, the common spatial 
frequency grid is associated with the equivalent full circular aperture SpinAp synthesized spatial frequency passband. 
[0050] Any or all possible distinct pairs of spatially registered and Fourier transformed individual frame images 30-4 
are multiplied by spatial frequency filters 31-4 designed to compensate for the temporally varying system transfer 

50 functions associated with each individual frame and resulting from the rotation of the SpinAp sensor system 10. The 
system transfer function compensation may be implemented by one of several methods, a few examples of which 
include: equalizing each frame to a target full circular aperture or other reference system transfer function, equalizing 
one frame to the system transfer function of the frame to be differenced, or deconvolving the system transfer function 
from both images 30-4 to be differenced. For ideal system transfer functions, and if and only if only frames acquired n 

55 radians apart are to be differenced, no system transfer function equalization is required. 

[0051] The equalized individual frames are differenced 40-4, and spatial frequency filtered 41-4 to eliminate non- 
common information outside the overlap passband associated with the two frames. The spatial frequency filtered frame 
differences are inverse Fourier transformed 42-4 producing individual pairwise differenced images 45-4. If the number 
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of distinct individual frame images is N, then the total number of possible distinct differenced individual frame image 
pairs can be as large as N(N-1 )/2. 

[0052] For a template signal detection approach, SpinAp moving object detection methods 20-4 based upon any 
one of the many possible equalization procedures have unique noise performance characterized by the particular 
5 method's receiver operating characteristics (ROC). The differences in ROC for signal template approaches are dom- 
inated by the different and uncompensated levels of noise "boosting" or "de-boosting" generated by the equalization 
process. However, an appropriately designed noise optimal matched filter approach compensates for the noise char- 
acteristics in the equalization process. 

[0053] To obtain the best possible detection performance, the set of differenced spatial frequency filtered transformed 
10 images 45-4 is transferred to a noise optimal matched filter detection processor 50-4. The matched filter detection 
processor 50-4 incorporates an array of individual pairwise differenced frame matched filters 55-4 having properties 
corresponding to the physical phenomenology of potential objects, such as, shape, size, velocity, acceleration, back- 
ground to object contrast ratio, signal to noise ratio, number of objects, number of track crossings, and the like. The 
matched filter detection processor 50-4 also incorporates a priori knowledge of the measurement noise statistics prop- 
15 agated through the equalization, differencing, and filtering process, as well as unique signal measurement character- 
istics of the SpinAp sensor system 10. Each of the individual pairwise differenced frame matched filters 55-4 of the 
detection processor 50-4 is applied to every pixel location for the appropriate frame difference and subsequently 
summed. The output of each object phenomenology total matched filter 56-4 is transfe'rred to a comparator 604 with 
thresholds set to a predetermined level for each of the total matched filters 56-4. A moving object detection is declared 
20 when any of the outputs from the total matched filter array 55-4 exceed the predetermined thresholds of the comparator 
60-4. The threshold may be predetermined by a priori assumptions concerning the signal and noise statistics, or by a 
posteriori assumptions concerning the signal and noise statistics based upon the individual frame or set of SpinAp 
frames. 

[0054] The equalization, differencing, and filtering operations eliminate all temporally stationary components in the 
25 common overlap passband of the individual frame images, leaving only the signature of smeared moving objects. The 
performance of the moving object detection method 20-4 is scenario dependent and is degraded by individual frame 
measurement noise propagated through the equalization, differencing, and filtering operations, and background leak- 
age due to residual frame to frame registration error. If observational parameters permit, acceptable performance may 
be achieved for less than optimal approaches, whereby the optimal array of matched filters 56-4 and comparator 60-4 
30 may be replaced by template matched filters and the comparator 60-4, streak detectors and the comparator 60-4, or 
a simple single pixel or multiple pixel threshold exceedence detection processor. 

[0055] For the purposes of completeness, the formulation of the SpinAp moving object detection processing method 
20 is as follows. Using the definitions and nomenclature in U.S. Patent No. 5,243,351 , a description of the mathematical 
foundations of the enhanced moving object detection SpinAp sensor 10 and data processing methods 20, 20-2, 20-3, 
35 20-4 embodied in the present invention is presented. 

[0056] Individual SpinAp noise free measurements will first be discussed. In the absence of noise, the mean number 
of detection electrons produced during the integration time of a detector located at position R acquired during frame 
time tj may be expressed as, 



40 



45 



50 



55 



Os(R,tj) = 1 d^OsjF.tjJexpt -2«if R I . 

where O s (F ,tj) is the spatial Fourier transform of the detector output at spatial frequency F and frame time tj. Assuming 
the effective sensor system transfer function is temporally stationary during the detectors integration time AT, the Fourier 
transform of the individual frames detected intensity may be expressed as 

I r*" \ 

Qfotj) .« lAsWoAfpasSTFtF.tj) U_ dtS n (MF,t) . 

I AT M / 

A 

where \ AS is the intensity spatially averaged over the effective focal plane array area, W 0 is the zero spatial frequency 
component of the detector transfer function, A^s is the effective focal plane array area, STF (F ,tj) is the total sensor 
system transfer function correspond to frame time tj, and S n (MF ,t) is the normalized spatial Fourier transform of the 
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geometrical projection of the scene at time t. 

[0057] A mathematical representation of a scene and the scene's geometrical projection to the imaging plane, in- 
cluding background, a single moving object, and scattering component will now be discussed. The intensity emitted 
from a scene containing a single object moving in front (relative to the SpinAp sensor) of a background, and behind 
transmissive scattering layers may be expressed as, 

S(r,t) = T(r,t) + B(r.t) - B(r,t)W t W T (r,t) + H(r,t), 

where T(r,t) is the intensity distribution of the moving object at scene position r and time t; B(r ,t) is the intensity distri- 
bution of the background at scene position r and time t, H(r,t) is the intensity distribution of transmissive scattering 
layers corresponding to scene position r and time t, and W T (r,t) is the object's inverse silhouette at scene position r 
and time t. The object's inverse silhouette is given by 

W T (r,t) = {1 if T(r.t) > 0 and 0 otherwise. 
[0058] Subtracting B(r ,t)W t (r,t) in the equation for S(r,t) corresponds to removing the background photons obscured 
by the moving object at scene position r and time t. 

[0059] Therefore, the geometrical projection of the scene's intensity to the focal plane of. the SpinAp optical system 
is represented by 

l g (R)=^S(R/M,t) = 
-^fT(R/M,t) + B(R/M,t) - B(R/M,t)W T (R/M,t) + H(R/M,t)], 



where M is the magnification of the optical system. Fourier transforming the geometrical projection and applying the 
Fourier transform similarity theorem yields, 

Ig(F) a FT(l g { R)] = FTj^Ls(R>M,t)) = S(MF.t) 

or 

-MFT{ltR/M,t)} + FT{B(R/M,t)} +.FT{H{R/M,t)) - FT{B{lVM.c)\VT{R/M,t))} , 



where FT and A refer to the Fourier transform operation. Application of the Fourier transform convolution theorem, and 
transforming variables permits the evaluation of FT{B(R/M,t)W T (R/M.t)}/M 2 yielding, 



^FT{B(R>M > t)W 1 {l^M f t)} = M 2 1 d^I^F-Fj.tjw^MF ,t) . 
M J 

[0060] The corresponding geometrical projection of the scene spectrum and its zero spatial frequency component 
are given by 

Ig(F) = f(MF,t) + H(MF,0 + B(MF,t) - M 2 J d 2 F'S(N{F-F ItJw^MF ,t) 
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and 



Ig(F = 0) a T(MF = 0,t) + H(MF = O.t) + B(MF = 0,t) - M*J d 2 F §(-MF '. t)W|(MF \t) . 
[0061] The instantaneous normalized geometrically projected scene spectrum is represented as 

Ig(F = 0) S(MF = 0) 
T(MF.t) + H(MF,t) + B(MF,t) - M 2 | d 2 FBMF-F],t)w^MF ,t) 

T(MF = 0,t) + H(MF = 0,t) + B(MF = 0,t) - M 2 J d 2 F §(-MF '.t^MF ,t) 

[0062] The normalized scene spectrum can be decomposed into two components dominated by the background and 
foreground spectrum, and the moving objects contribution to the spectrum, such a representation is 



S n (MF,t) m B n (MF,t) + AT„(MF.t) 



where 



B„(MF,t) = 
fi(MF,t)+§(MF,t) 



[ft(MF = 0,t) + B(MF = O.t)] \ I + - 



T(MF = 0,t) - M 2 j d 2 F §(-F M.tJw-^F M,t) 



H(MF = O.t) + B(MF = O.t) 

A _ 

[0063] B n (MF ,t) is the normalized background and foreground component for spatial frequency MF at time t, and 



T{MF,t) ^ M 2 j d 2 F'B([F-FlM,t)w 1 {FM f t) 



) T(MF = at) - M 2 1 d 2 F§(-F Mu)Wt{f M,t) 
[fi(MF = 0,t) + B(MF = 0,t)]\ 1 + ——J——— . 



H(MF = 0,t) + B(MF = 0.t) 



where AT n (MF ,t) is the normalized object's variation from the background and foreground for spatial frequency F at 
time t. 

[0064] Normalized background and foreground spectra for a temporally stationary scenes will now be discussed. 
For observational scenarios having temporally stationary total background components, the normalized total back- 
ground spectrum, and moving object variation are 
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B^MF.t)s _ 
^MF,t) + B(MF,t) 



(S(MF = 0,t) + B(MF = 0,t)] \ 1 + • 



T(MF = 0.t)-M 2 



d 2 F'S(-FM,t)WT(FM.t) 



H(MF = O.t) + B{MF = 0,t) 



and 



AT„{MF,t)s 
T(MF.t) - d 2 F'S([M(F.F )])Wf(F'M,t) 



(H(MF = 0) + B(MF = O)] 



1 + 



T(MF = O.t) - M 2 J d 2 F'S(-F 'm)wJf M.t) 
H(MF = 0) + b(mf = 0) 



[0065] Moving object and background decoupling will now be discussed. If the number of background photons ob- 
scured by the moving object is negligible when compared to the total number of background photons, then the total 
background and moving object variation normalized spectra can be further simplified, yielding 



B n (MF,t) ^ B n (MF) - i B(MF)+ § (MF) 
" " H(MF = 0) + B(MF = 



0) 



and AT n is given by 



~ _ T(MF;t) - M 2 1 d 2 F'§([MF-F ])wj(F M.t) 



a¥„(mf.0 = 



[h(mf=o) + b(mf=o)] 



[0066] Temporal variation of the object and the object's inverse silhouette will now be discussed. The time dependent 
motion of the object may be describe by applying a temporal dependent spatial shift to the argument of the objects 
spatial intensity distribution. Defining the starting time of the imaging process to be to, the initial object and the object 
inverse silhouette are given by 

T 0 (r) = T(F,t 0 ) 



W T0 (r) = W T ("r,t 0 ). 

[0067] The motion dependent temporal variations of the object and the object silhouette are given by 

T("r.t) = T 0 fr-r T (t)) 
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W T (?,t) = W T0 (P-? T (t)). 



T(MF.O = ^MF,to)exp(.2«iMF FKO) 
a T<j(lvff)exp(-2itiF R-Kt)) 



10 



15 



and 



Wt{MF,i) = WT{MF,to)exp(-2niMFri(t)) 
9 W-rt(N#)exp(-27tiFRi<t)) 

20 

25 yields 

at„([m(f-f)])= 

Td(Nff)exp (-2*iFR T W) - M*J d^B( Mf-F )] )WTo(FM)exp(-2^F RT(0) 

[h(mf = o)+b(n#=o)] 



55 =^^^r^r^=^^.^^^ b = 

variation, AT n , 

40 AT n (MF,t) = 

T(j(Nff)exp(-2rt iF V(t)) - M^J d^FB(^F-F])w^Nff W-2«g V(t)) 
[fi(Nff = 0) + B(lvff = 0)] 



50 



55 
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10 



15 



20 



45 



50 



OstF.tj) = iAsWoAfpasSTFfF.tj) (b n (MF) + ->-f diAT n (MF.t)l 

= iAsWoAfpasSTlKF.tjjBMlNff) + lAsWoAfpasSTFfctj) • 

To(MF) sinUF-VATj^ .^p,^,) 

[H(Nff = 0) + B(MF = 0)] TiFVAT 
-iAsWoAfpasSTFtF.tj) • _ 
M2 [ , F b(m[f-f])Wti/mf) sinfoFVAT)^ .^ 
J [fi(MF = 0) + B(Nff = 0)] TiFVAT 

where the time average of the constant velocity motion required integration of the 'track' integral 



f * ._. exn(-2itiF- V t) , 
- dtexp23tiF VT = -i2j —I I _^ 

J |jf cxp(-27tiFV) *¥ 

, \sin (hFVAT) 
s exp(2TriF Vtj) _ _ -. 



25 nFVAT 

[0072] The SpinAp synthesizer and synthesizing process will now be discussed. The simplest version of the SpinAp 
synthesizing process employs a synthesizer, such that the spatial frequency redundancy average of the noise free 
30 measurements yield the equivalent full circular aperture image, or 

0 Fsyn (F)=A<F)6 Save <F) 

A _ 

35 where 6 F syn (F ) is the full circular aperture image to be synthesized, A(F ) is the synthesizer, and where O s ave (F ) is 
the spatial frequency redundancy average of the (noise free) individual frame measurement spectra. 

6 Savc (F) = -LrX 6^F,tj). 



[0073] The synthesizer, A, is given by 



^ iAFAfpaFWoFW^F^TF^F) _ WAfoapWoFSTF^ 
UiAfpasWos W NS OTFs(F) UsAfpasWosSTF^F)' 



f 0r W NS OTF s (F) = STF S (F) * 0, where STF corresponds to the system transfer /unction, which is the product of the 
optical system transfer function, OTF, the normalized detector transfer function W N , and the subscripts S and F des- 
ignate parameters associated with the SpinAp strip aperture sensor and the target full circular aperture system, re- 
55 spectively. The spatial frequency redundancy averaged strip aperture transfer function appearing in the denominator 
of the synthesizer is 
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W NS OTFs(F) = STFs(F) = 

Nrf(Fjj = i j=l N rf (FJ 



[0074] Therefore, the synthesized image spectaim resulting from the observation of a scene consisting of a tempo- 
rally stationary total background and foreground, and an object moving in the sensors field of view is given by 



n^f) 

n,.(f) _ ^ 

= A(F)— ^ X lAsWoA fpa sSTF(F,t j )B N (lvff) + 

N rj (F) j = | 

+ A (F>-J_T lAsWoA^asSTFtF.tj )^^ _ _ 

H<(F)j = i h(mf = o) + b(mf = o) 

sin(«F^VAT) exp( , 27t . F . v ^ 
JtFVAT 

-A(F}— L=r N X lAsWoAfpasSTFtF.tjjM 2 • 
N r t(F) j= , 

d2p . §(m[f-F ])Wtc(mF ') sinfoF ^T)^ ^ . y ^ 
H(MF = 0)+B(MF = 0) tiF VjAT 



[0075] Further simplification produces, 



Of syri(F) = iAFWoFAfpaFSTF^FlB^Nff) + 

j. T w a STF^F] Tp(MF) 

+ lAFW ° FAfpai ra fi(MF = 0) + B(MF = 0) 

£ STFstF.ti) 

j.l TtF-VjAT 

lAF W 0FAfpaFSTFf{F) N ^ STFs(F,t ) • 

I STFs(F,tj) 
j = » 

M 2 [ d y BMF-FpW^Nff ) sin(rf jVjAT)^ ^ • ^ 
I H(MF = 0) + fi(MF = 0) tiF VjAT 
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[0076] The first term in the synthesized image spectrum equation corresponds to the equivalent full circular aperture 
image of the temporally stationary background. The next two terms in the equation correspond to the moving object's 
signature in the synthesized image. 

[0077] Differenced synthesized images and moving object signatures will now be discussed. A moving object spatial 
frequency domain signature may be obtained by differencing two or more synthesized image spectra. Differencing two 
completely independent, full aperture synthesized images (or image spectra), while observing the same scene, requires 
a minimum of two 180 degree rotations. Alternatively, to decrease the time between frame differences, synthesized 
images (or spectra) may be differenced with as few as one non-common frame, or the synthesis operation may be 
performed with fewer than the required number of frames to generate the full circular passband image spectra. 
[0078] Since the temporally stationary parts of the synthesized images are identical for any pair of synthesized images 
observing the same scene, the difference of synthesized images will contain only moving object information. The noise 
free spatial frequency domain moving object signature of the differenced synthesized images is given by 

AOba(F) '- 6 jyn „(F) - 6^ A (F) 
lAFWoFAfpaFSTFFtF) Tq(MF) 

NttfO _ H(MF = 0) + §(lvff = o) 
£ STF s (F,tj) 

^ jTi JlFVjAT 

*M ( d y §MF-f])Wt^MF) . 

- £ STFs(F,tk)M 2 | H(MF = 0) + B(MF = 0) 
J 

TCF -V k AT 

[0079] The spatial signature of the moving object resulting from the difference of synthesized images is 




a6ba(R,St)= I d 2 FAOBA(F)exp(2jriFR), 



where 5T represents the time delay associated with the acquisition of different sets of acquired strip aperture frames. 
[0080] Differencing any pair of appropriately modified SpinAp frames for background removal and moving object 
detection will now be discussed. For 2Nf individual SpinAp strip aperture frames, there are 2Nf (2N f - 1)/2 possible 
distinct pairwise frame comparisons, or differences. An optimal detection technique would take advantage of the back- 
ground removal when differencing any possible pair, as well as, the potentially large number of moving object pairwise 
signatures. However, unlike differencing frames for the purposes of background removal and moving object detection 
for temporally stationary aperture systems, specific techniques and methods must be applied to the SpinAp data sets 
to compensate for the rotating nature of the system transfer function. The SpinAp data sets are modified to distinguish 
between temporal variations in the scene and temporal variations originating from the moving object. 
[0081] SpinAp individual noisy frame measurements and frame to frame system transfer function equalization will 
now be discussed. The individual strip aperture noisy measurements spectra can be represented as the noise free 
strip aperture detected image spectra plus the strip aperture additive noise spectra 
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C S m(f : -*i) = 6 s(F.V + "s( F .«i)- 

A 

where C s and n s are the noise free signal and noise spectra, respectively. 

5 [0082] Individual frame system transfer function passband filter will now be discussed. In general, the individual 
frame measurements provide information related to different portions of the equivalent full circular aperture optical 
system passband. The detection of temporally nonstationary phenomenology, such as, moving objects or transient 
events requires the identification and utilization of common information between frames, as well as, accommodating 
for the differences in image contrast due to the time varying nature of the SpinAp sensor's transfer function. The frame 

10 to frame system transfer function compensation can be implemented by several methods, a few examples of which 
are equalizing each frame to a target full circular aperture or other reference system transfer function; equalizing one 
frame to the system transfer function of the frame to be differenced; deconvolving the system transfer function from 
both images to be differenced; or for ideal system transfer functions, and if only frames acquired n radians apart are 
to be differenced, no system transfer function equalization is required. 

15 [0083] For a template signal detection approach the SpinAp moving object detection methods based upon any one 
of the many possible equalization procedures would have unique noise performance characterized by the particular 
method's receiver operating characteristics (ROC). The differences in ROC for signal template approaches are dom- 
inated by the different and uncompensated levels of noise "boosting" or "de-boosting" generated by the equalization 
process. However, an appropriately designed noise optimal matched filter approach would compensate for the noise 

20 characteristics in the equalization process. 

[0084] While any one of many methods can be used to implement the system transfer function equalization, the 
formulation described in the following paragraphs is typical of the mathematical details associated with each procedure. 
For brevity, only the procedure whereby one frame is equalized to another is described in detail. 
[0085] Selecting one frame as the current reference frame and any different second -frame to be used in the com- 

25 parison, the spectra of the frames are spatial frequency filtered to provide spatial frequency information in the common 
overlap region of the passbands of both frames. The second frame is also modified to compensate for the differences 
in system transfer function, i.e. 



30 



m R (F,ti) = xtF.tjMF.tipsmfF.ti) = xfF.tj&F.tjforfF.tj) + n s (F.tj)], 
and 

X(F,ti£(F, tjlOsiF.tj) + ^(F.tjjjHtF.tj, 

A A 

where m R (F,tj) denotes the spatial frequency filtered output of the reference Jrame, and m E (F,tj) the spatial frequency 
to filtered and transfer function equalized output of the frame to be compared, x is the optical system passband filter, 

J(F,t j )=1ifSTF s (F,t i )>0, 

45 and Q otherwise 

and H is the noise free passband equalization filter 



50 



55 



STF s (F t tj) ^ 
H ( F 'V^= 1lf STr^ >0 - 

and 0 otherwise. 

[0086] The filtered and equalized difference between any two frames can be expressed as the sum of differenced 
signal and noise spectra, i.e. 



Amij{F) « m R (F.ti) - rr^F.tj) = AS./F) + AN|/F), 
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where the differenced, filtered, and equalized noise spectrum, an m , is 

ANj/F) = x(F.tjWF.ti)ns(F.ti) - xtF.tiMF.tjjnsfF.tjJHsfF.tj.tiX 

and the differenced, filtered, and equalized signal spectrum, 5* (F), is 

ASij(F) = lAsWosAfp a sX(F,ti]X{F,tj)STF s (F,tis N {M,F,ti)] 
= I^WosA fpa sXjF.ti)X(F.t J )STF s (F.ti) • _ 
• [§N(M.F,ti) + T n(M,F,ti) - BNtM.F.tj) - T^M.F.tj)]. 

[0087] For temporally stationary backgrounds and foregrounds the resultant/signal spectrum is related to the differ- 
ence between the normalized moving object variation from the background, AT N , by 

^ij|F) = T A sW 0S Afp a sX(F.t i WF,t j )STF s {F,ti(T N {M,F,t i ) - AT^M.F.tj)].. 

[0088] If the object is moving with constant velocity during the integration time of the detectors, then the filtered 
differenced signal corresponds to 

ASij(F) - lAsWosAfpasxtF.tiJxtF.tjjSTFsfF.ti) • 

%ML- sin(nF lV T) [exp(.2 n iF. Vfo) - exp(-2niF- %| - 

H(MF = 0) + B(MF - 0) «F • VjAT 

+ lAsWosAfpasxfF.ti^CF.tjJSTFsfF.ti) • 
I I H(MF = 0) + B(MF = 0) tiF'-VAT 

J H(lvff = 0) + B(Nff = 0) nF-VAT J 

[0089] SpinAp spatial frequency and temporal frequency domain moving object detection will now be discussed. A 
simple mathematical representation of the spatial frequency and temporal frequency domain SpinAp moving object 
detection procedure is obtained by approximating the discrete temporal sampling as a continuum and taking the con- 
tinuous Fourier transform of the spatial frequency transform for each of the SpinAp individual frames. The Fourier 
spectrum of the noiseless individual frame measurements is given by 
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OstF.tj) = TAsWoAfpasSTFsJF.tiJBMMJ) 

P H(MF = 0)+B(MF = 0) nFVAT. 

.[AsWoAfpasSTFtF.tj). 

I H(Nff = 0) + B(MF - 0) irFVAT 



[0090] The temporal Fourier transform of the spatial frequency spectrum of the individual frame is obtained by eval- 
uating the integration, and applying the convolution theorem yielding 



S(F,f) = | dt exp(-2nift)6s(F,t) 
= UsWoArpas j dt exp(-2nift)STFs(F,t)S > j(lvff ,t) 
= lAsWoA rpa sj df STF s (F,f)SN(lvff,f-d 



where Stf s is the temporal Fourier transform of the SpinAp system transfer function, anda is the temporal and spatial 
frequency Fourier transforms of the normalized geometrically projected scene spectrum. 

[0091] Evaluating the temporal Fourier transform of the normalized geometrically projected scene spectrum, the 
noiseless temporal frequency transform of the strip aperture measurements, becomes 

6s(F,f) = I A sWoAfpaS ' 

+ A «* f pjV • jgw WM.F) sin fe 1 - AT) STF s (F,f-F- V ) + 

AfpalHfpa + Bfpa) JlF-VAT 
Aug fpaBfpa 



Afpa(Hf pa + Bfpa) 
dywWMl f ( ?'5 AT) A B N d(N{F-F ])STFs(F,f-F V ) 



nF VAT 



[0092] The first equation above defining C s (F,tj) is used to determine the signal component of the temporal and 
spatial frequency matched filters. The matched filter essentially eliminates the temporally stationary components, since 
the first term of equation for ^(F.f) contains all temporally stationary effects and is zero for all nonzero spatial frequen- 
cies. 

[0093] Simplified signature's for uniform background, foreground, and moving object will now be discussed. Predic- 
tion of receiver operating characteristics for SpinAp moving object detection procedures can be simplified for the special 
case of observational scenarios associated with uniform intensity backgrounds, foregrounds, and moving objects. In- 
deed, the design of coarse moving object detection approaches may utilize such simplifying assumptions. The tem- 
porally stationary and spatially uniform components of the background and foreground can be expressed as. 
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B(r.t) = B 0 s B SCN and H(r,t) = H 0 ■ H SCN 

■ where B SCN is the average background intensity in scene coordinates, and H SCN is the average foreground intensity 
in scene coordinates. * 
[0094] Focal plane average values and Fourier transforms will now be discussed. The geometncal projection of tne 
average values of the scene's background and foreground intensities are and 



and the corresponding Fourier transforms are 

B(MF) = B^8(MF) = ^B^S(F) = B^S(F) 
M 



H(MF) - H scn 8(MF)= ^H scn 5(F)= H^F) 

where use has been made of the Fourier transforms of the backgrounds in the scene coordinate system 

B(F) - B" 0 8(F) = B^5(F), and H(P.t) - H~ 0 8(F) = H^8(F) 
and the Fourier transform similarity theorem. 

10095] Zero spatial frequency relationships to mean values will now be discussed. The normalized geometncal pro- 
jection of the scene spectrum is generated by dividing the Fourier transform of the geometrical projection of the scene 
by the zero spatial frequency component of the geometrically projected scene spectrum. The geometrically projected 
zero spatial frequency component of the background's, and moving object's intensities can be expressed in terms of 
the spatially averaged mean values and the geometrical projected area of the backgrounds and object, as 

T(MF=0,t) = TfpaA^^, B(MF=0 ( t) = B^A^. H(MF=0,t) = H^A^, 

where At rqfpa is the approximate instantaneous geometrically projected area of the moving object, is the approx- 
imate detected area of the background in the focal plane of the SpinAp sensor, and the overiine symbol (-) refers to a 
spatial average in the focal plane array. 

[0096] Zero spatial frequency relationships for the object inverse silhouette will now be discussed. Applying the 
Fourier transform definition to the object inverse silhouette in the scene coordinates, and transforming to the image 
plane coordinates by a change of variables given by Mr scene = R lma g e . yields 



Wt(MF,i) = I d 2 rW-KM)exp(-27tiMFr) = 

/seem 

-J d 2 RWj(R/M,t)exp{-27CiFR). 



M 2 J 

[0097] Therefore, the geometrical projection of the moving object's zero spatial frequency component can also be 
related to the object's effective area in the focal plane by 
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. .-WT<MF=o,t.)-^ j-d?Rw-KRM,o = ^aja. - 

5 

[0098] Normalized spectra for constant foreground and background will now be discussed. Applying the preceding 
relationships for the uniform background components of the scene enables determination of mathematical expression 
for the normalized background and normalized object. Recognizing the ratio of the sum of background and foreground 
10 projected scene spectra to their zero spatial frequency component to be 



H(MF) + B(MF) _ H fp , + B fpa S(F) _ 5(F) 
H(MF=0) + B(MF=0) ~ (Hi^ + SfSJAflpa Af P a 



30 



yields the total normalized scene background given 

iMMF,t) s §n(mf) ■ ir~i — " — ~k — =u- 

\ Afpa(iV.+ Bl^)/ 

A 

25 [0099] In like manner, the normalized object difference from the background, T N , for temporally stationary and spa- 
tially uniform backgrounds can be expressed as 

~ x A fp a (H^ + B^J J A tpjHfpa + Bf pa ) 

ATn(m,f,)s 

Af pa (Hfp a + Bf pa ) 

[0100] Uniform object approximation will now be discussed. In scene coordinates, a uniform object can be repre- 
sented as a product of the object's inverse silhouette, W T , with the object's constant intensity level, T 0 , i.e. 

T(r,t) = T 0 W T (r,t). 
Therefore, the Fourier transform of a uniform object is 

T(F,t) = T scn W T (F,t), 
and the corresponding geometrically projected transform of the uniform object is 

T(M,F,t) = T scn W T (MF,t). 

A 

[0101] Thus, for uniform objects and backgrounds, the normalized object variation from the background, AT n , for 
temporally stationary scenes is given by 
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TijMF.Q _ Ml f d 2 F - B^^ F - F JL Wt{ivF \t) 
AT„(MF,t) = - JSLJ 

A fpa(Hf pa + Bfpa) 



_ (Tj^ - B^) Au, fpa WiN(MF,t) 

\ Ajp^H^-i-B^ / 



where the normalized inverse silhouette is defined to be 

WiM(MF=0,t) ±S1 55 
M 2 

and the zero spatial frequency component of the normalized inverse silhouette is 

W T (MF=0,t)=^2|2. 

M 

[0102] Applying the Fourier transform shift theorem to describe an object moving at constant velocity the object 
inverse silhouette is 

W TN (MF,f) = W TN (MF,t 0 )e X p(-2niFVt) = W TNo (M,F)exp(-2niFVt) 

where W TNo (MF ) = W-^MF .y is the normalized object's inverse silhouette's Fourie transform at time = ^ The resulting 
expression for the object variation from the backgrounds is 

A * | M ? ( )- ^ - ^ A ' rg fpaW WM.FVxpj^TCiF^Vt) 
AT„(M,F,t)- ( __- +5; p Afpa K t Atrg (pafTfpa ~ Bfpa)\ 

\ Afpa(?V + B^ / 

[0103] The differenced and filtered moving object spatial frequency signature for uniform objects and back- 
grounds^, is therefore 

ASiiflF) S I A sWo S Afp a sX(F,tiMF.t J )STFs(F.t i (saMF,tO^MF.t j )] 
~ r~ lr - s^.n >sin( : ffF V'AT)(T^-B^) 

j , A tf g f P a(Tfpa - Bfpa) 
Afp^Hfpa + Bfpa) 
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[01 04] The spatial domain differenced signal for uniform object, background, and foreground will now be discussed. 
Inverse Fourier transforming the moving object's differenced and filtered spatial frequency spectrum, produces the 
spatial domain moving object signature, which is 



-ASij(R) = j d 2 F exp(-2*iF. R }*S { 0) 



_ ^ -FasWqs Atrg fpa (Tfoa * Bf pa J 
j + Atrgfp^TfpaS " Bf pa s)( Hfpa + Bfpa) 
Afpa(HfpaS + Bfpas) 

f m , * ,sin (tcF- V AT) fjy - B^) 

J d F exp(-2niFR K(FAte(F,tj)STFs(F.ti) ^ ( gg^ - 

[01 05] The differenced detector mean signal to RMS noise ratio will now be discussed. The mean differenced signal 
to RMS differenced noise ratio at location R for uniform objects and backgrounds can be expressed as 



25 



30 



35 



40 



45 



50 



SNRij(R) = 



ASjj(R m ) _AS ij (R m ) 



7 u> a (TfpaS - Bfpas) 

VA sS _22__2s2_^r/i , ATOfpalTfpas - Bfpas A 
4 D x Ax s D y Ay s I AfpJHfpaS -i-Bfpas) / 

| ^Fexpl^RiF R^F.ti^F.tjJSTFjKF.ti) • 
sih(^VAT)-^ I ^ e . 2jliF . v(ti) e . 2KiF v(,,)} 



TtFVAT 



V | d 2 F exp(-2niF- (R m - SJ) ( 1 + |8(F,tj,tjj 2 jriF.tj&F.tj) 

where <-> denotes noise ensemble averaging, and the differenced noise variance oJ^RJ ma Y be obtained by 
analytic procedures described with reference to spatially white noise. 

[0106] Any differenced synthesized pair signal spectrum for uniform objects and backgrounds is now represented. 
The noise free spatial frequency domain signature for any pair of differenced synthesized images associated with 
uniform objects and backgrounds is given by 



55 
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- t-\ 7Z ic\ ^WqfA fp f(F) Atrs fpaCTfpaS - B fpa s) 
O iynBA (F) = O syn8 (F)- OsynAfF)^ aJh^I + B^) 

X STFs(F,tji) 
j = » 

\ k= i nF -V k AT / 

and the corresponding spatial domain moving object of any two differenced synthesized images is 

AOba(R.5T) = j d 2 F OB A (F>xp(-2jiiF R ) 

AF OF fjpF Afpa(Hf p aS+Bfpas) 
• ( d^Fexp^iFR ^^^ WtooCmf) 
^ X STFs^.tji) 

. STT^t r ^l^e xpf^iF-V^,)) 
^ jcFVjAT 

. N f STF 5 (F,t k r^3 AT L p(.2mF.V k (t.)) 
k=i wFVicAT 

r «> a A t f g fpjTfpaS - Bfp»s) . 
= lAsWosA^ fpaAf ^ j _ + g -_ ) 

[ j2c _/ • .= ^- vsin(jcF V k AT) STFf(f) 

I d 2 Fexp(2«iF R r ~T~"7?v 

J nF V k AT 

' X STFs (F.tj) 

i-i 

|N,r(F) \ 

• | X STF s (F.t>xp(.2JtiF.y]j(t J )) |w TO0 (Nff) + 

f N„(f) K 

- £ STT^F.tfJcxp(-2JCiF-yi(t k j WTOo(Nff). 

[01071 Matched filter nomenclature will now be discussed. The best possible performing detection processor, for 
wide classes of performance criteria, has been established to be an array of generalized matched filters with a com- 
parison device, which will be referred to as a comparator. The design of the optimum matched filter arrays for Sp.nAp 
moving object detection processors requires incorporating knowledge of the unique signal measurement charactenstics 
of the SpinAp sensor, the desired data processing before input into the matched filter array, and the an .opated object 
phenomenology. The a priori object knowledge can be either statistical or deterministic in nature, .nclud.ng physical 
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parameters, such as, shape, size, velocity, acceleration, background to object contrast ratio, signal to noise ratio, 
number of objects, number of track crossings, noise cross-covariance, etc.. The matched filters of the array for a 
particular signature can be expressed in matrix notation as, 

5 F - S*C- 1 

where S is the a priori noiseless signal vector, and C^ N is the noise cross-covariance matrix inverse. Both the signal 
and noise cross-covariance matrices are appropriately defined for the desired moving object detection data process.ng 

[uTo8] d The signal matrix resulting from any of the combinations of differenced frames, differenced synthesized im- 
ages, temporal frequency variations or other signatures as determined to be appropriate for SpinAp mov.ng object 
detection may be represented as 



15 



20 



\ s 



Si 



25 



where | refers to any or all combinations of independent spatial, spatial frequency, temporal, and temporal frequency 
signal coordinates associated with the a priori statistical or deterministic SpinAp moving object signature. 
[01 09] Likewise, the noise matrix resulting from any of the combinations of differenced frames, differenced synthe- 
sized images, temporal frequency variations or other appropriate signatures as determined to be appropnate for SpinAp 
moving object detection can be represented as 



30 



35 



N = 



Ni 
N 2 



N 



nRi) 



/ 



[0110] The set of noisy measurements can therefore be represented in matrix notation as the sum of signal and 
noise components, yielding 



50 



55 



M = 



Mi 
M 2 











+ 











[0111] Complex functioned matched filters can be constructed by recognizing the signal and no.se matrix elements 
are the real and imaginary parts (or amplitude and phase parts) of the complex signal and noise, respectively. 
[0112] For example, the signal and noise matrices for any number of differenced frames, or differenced synthesized 
images within the spatially shifting matched filter acquisition "window" would be represented by 
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/ 



s « \ 



'^iignil 



AS 12 (Ri) 

AS, 3 (R no J 
AS oP (R,) 
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10 



15 



20 



25 



35 



Ni 



N = 



AS I3 (R,) 
AN 13 (R2) 

AN a p(Ri) 
AN a fj(R*J / 

30 respectively, where n now refers to the number of total elements in the spatial matched filter acquisition window. 

[0113] The matched filter output: The matched filter output therefore corresponds to the noise free matched filter 
signal and additive "whitened" noise, or 



AN l2 {R|) 
ANu(R 2 ) 



Mout = FM = stC NN< S + N )- 



[0114] Matched filter receiver operating characteristics for Gaussian noise will now be discussed. The receiver op- 
erating characteristics (ROC) for the SpinAp matched filter moving object detection processing methods utilizing a 
comparator and a threshold exceedence decision approach can be readily determined for approximately Gaussian 
individual frame measurement noise. The probability of detection and probability of false alarms for the SpinAp matched 
filter output for a given threshold can be determined from 
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2 ^ Vis^ci; 

where In K is the threshold parameter of the comparator, and 



C WL+... InJK. 



2fl ~ V2S , Cn' n S 



) 



respectively. 

101 1 51 SpinAp differenced noise cross-covariance matrix elements will now be discussed. Determination of the noise 
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optimal SpinAp moving object detection matched filters depends upon a priori knowledge or assumptions concerning 
the noise processes and noise propagation from the measurement steps to the matched filter. The classes of SpinAp 
moving object detection procedures described elsewhere in this disclosure utilize individual weighted frame differences, 
synthesized frame differences, or temporal and spatial frequency matched filters. The determination of the appropriate 
propagated noise cross-covariance matrix elements is required to establish the noise optimal matched filters, or the 
performance of less than noise optimal techniques. The following sections provide example techniques for determining 
the appropriate SpinAp noise covariance matrix elements. Covariance matrix elements associated with processing 
steps not explicitly described can follow similar methods and procedures of determination. 

[01 1 6] The noise spectrum for any pairwise differenced frames will now be discussed. One embodiment of the SpinAp 
moving object detection processing method differences any distinct pair of system transfer function equalized, and 
overlap passband filtered individual frames. The corresponding differenced, equalized, and filtered noise spectrum, 



Nij, is 



ANij(F) a JtF.tiMF.tiJnstF.ti) - xtF.ti^F.tjJnstF.tj^F.t^ti). 

[0117] The individual pixel noise associated with position R is obtained by performing an inverse Fourier transform 
of the processed noise spectrum, and is 



ANi/R) s j d 2 FANij(F)exp(27CiF.R). 
= | d 2 Fexp(27uF.R ^(F^ifctj^nstF^O * ^(F^ij^F,^,!,)} . 

[0118] For zero mean individual frame noise, 

(nsGU,)) = 0. (pM-hlj = 0 (nslF.ti)) = 0. 

(SstF.tj)) = 0. and (aN^R)} = 0 = (aN^F)) . 

[0119] The spatial noise cross-covariance will now be discussed. The spatial noise cross-covariance between any 
two pairs of differenced, equalized, and filtered frames at pixel locations R m and R n is 

(Ari/R m )Arki(R n )) = 
|d 2 Fjd 2 F'exp[-2jti(F R m -FR„)] • 

/H(F.t j ,t i )H>',t 1 ,t k )fis{F,t j Ks(F 4 + ^s(F.tiK{F .t fc )> - \ 
' \ - H*(F,t 1 ,t k )^s(F,t j )S' s (F,t 1 }- ^(F.tj.t^F.tjjn'siF'.tk)) / ' 
• x(F.tiK(F.tjK(F,t k K(F.t,) 

where the super script * refers to complex conjugation. 

[01 20] The effects of zero mean, temporally stationary and uncorrelated noise will now be discussed. For zero mean, 
temporally stationary, and temporally uncorrected noise, the individual frame noise spectrum cross covariance is re- 
lated to the noise power spectral density C nn by 
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(SsfF.tj^F-.ti)) = S i4 C ri ri(F)8(F - F), 

5 where 6j j and 5{F - F) are the Kronecker and Dirac delta functions, respectively. 

[0121] Substituting individual frame noise spectrum cross covariance into the expression for the spatial differenced 
and filtered noise covariance produces 



10 



15 



(AN ij(R m )£Nki{Rn)) >- j <1 2 fJ d 2 Fe-2p«(F Rn.-F' R„) . 

{H(F.ti,tj)H'(F,t,,t k ) 8 jt i + 8 jlk - H'fF.ti.tj^u - HCF.tj.tiJSj.k } 
• Zp.tiS(KtjK(F.t k K(F,t,)c ilI »(F)5(F-F) 



[0122] A special case corresponding to the pixel to pixel cross-covariance for a single filtered and equalized frame 
20 difference is obtained by letting the frame indices take the values i = k, j = 1 , and i ^ j. The corresponding effective 
noise spatial cross covariance is: 

(^Ny (R ro )ANij (R„)) = 
j d 2 Fexp(-2TuF- (R m - R n ))j 1 + 1 HfF.ti.tj)! 'feF.ti^.tjjCnntF). 

30 [0123] Spatially white noise and Nyquist sampling will now be discussed. For spatially white noise, and Nyquist 
sampling of the optics passband the individual SpinAp frame noise power spectral density can be approximated as 

A - _ A dS XZ XZ — 
C nn( F ) = 4Dx Ax s D y Ay s a ^ 

where X is the wavelength, Z is the effective optical system focal length, D x and D y are the aperture sizes in the x and 
y directions, respectively, Ax s and Ay S are spatial sampling increments, which for 100 % fill factor focal plane detector 
arrays corresponds to the detector spacing, and o£ s is the mean noise variance on an individual detector. Substituting 
40 individual white noise power spectral density cross covariance into the expression for the spatial differenced and filtered 
noise covariance produces 

(ANij(R ro )AN kl ^ n ))= *f *Z to d2Fe-*«F fo-R.) . 

* DxAxsDyAys ; 

• (H(F,t j .tj)H*(F,t 1 ,t w ) 5 j4 + 5i, k - H*(F.t,.t k ) 8 M - H(F,t,,tj) 5 iJt ) • 
5c(F,tiS(F,t j H(F,t k g(F,t l ). 
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[01 24] The special case corresponding to the pixel to pixel cross-covariance for a single filtered and equalized frame 
difference is 
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(AN ij (R n ,)AN;^R 0 )) 3 
^3_XZ__XZ_^[ d *Fexp(-2*iF • (R m -Rn))(l -| H^.t^ftF.^.g. 

[0125] Pi-wise differences and the covariance matrix elements and white noise will now be discussed. One further 
simplifying example, which for ideal system transfer functions has the processing advantage of not requiring system 
10 transfer function equalization, corresponds to differencing frames that are n radians apart. For an ideal rectangular 
strip aperture, the system transfer functions for two frames n radians apart frames provide identical contrast perform- 
ance, and therefore the equalization filter is exactly unity, or HfF.tj.tj) = 1, and the overlap passbands are identical. 
Therefore for spatially white, and temporally uncorrelated noise, the cross covariance noise matrix elements are 
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15 

(ANijfRmjAN^ilRn^CNNlij = 2<46 M 
or the covariance matrix can be represented by 

A _ — 

C NN - 2 °dS j 

where I is the identity matrix. The inverse of the cross covariance matrix is therefore 

^NN == 1 

2<4 

30 [0126] The synthesized noise spectrum will now be discussed. Another embodiment of the SpinAp moving object 
detection processing method differences any distinct pair of SpinAp synthesized images. As described in U.S. Patent 
No. 5,243.351 , synthesized noise spectrum for the simplest form of synthesis process is 

35 JJ»yn(F)=A(F)n Sa v e (F) 

A 

where the spatial frequency redundancy averaged noise n Save (F), is given by 

"Save<H=j^2: n s (F, fj). 



40 



[0127] The differenced synthesized noise spectrum will now be discussed. Given the synthesized noise spectrum, 
45 the noise spectrum, An syn (F), for any pair (not necessarily derived from completely distinct data sets) of synthesized 
images is 
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An syn (F) = AAnsav e (F). 

A 

where An Save (F) is 

^-""idpiS Ssa(f -' ,) 
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where the subscripts A and B refer to the SpinAp individual frame sets that were used to generate the synthesized 
images. 

[0128] Spatial. domain differenced synthesized noise will now be discussed. Performing the inverse Fourier transform 
of the differenced synthesized noise spectrum, generates the spatial domain synthesized noise at pixel location R , 

n sy „(R) = | d2Fn sy ^F)exp(2KiF R ) = J d 2 Fn s ave (F)A(F)exp(2JtiF • R ). 

[0129] Spatial domain differenced synthesized noise covariance will now be discussed. The pixel to pixel cross- 
covariance for the synthesized noise difference can be expressed in terms of the SpinAp synthesizer and the cross 
covariance of the redundancy averaged differenced noise spectrum, or 

(ANs y j(Rm)ANt vn (R n )) = 
^| d 2 FA(F)nsav e (!W2niF *m)j d 2 FA (F & ave (F )exp(-2jriF ' R n )^ 

d 2 F'exp 2ni¥ R ra - F' R n )A (F)fe »ve(F)n£ *Jf } 



-H 



where the redundancy averaged differenced noise spectrum is 



rrWZ as B (F,t J )--J_5: Hs A (F,t k ) 
N rf (F) j=l N rf lF) k= , 

^Sav^/nSavdf// = \ ( ^) Ni /p) 



N rf (F)it; N rf (F) k =. ) 



nUf) n„(f") 

— N^FlNriCF)^!-!^! 

n^f)i4f), 

+ — -f-prl Z \Ssa(F, t k )nS A (F.l m ]} + 
n^f) hJf) 

NrftFjN^F)^! 1=1 

- * H I Z ^SA(F,t k )rrsB(F,a. 

[0130] Zero mean, temporally stationary and uncorrected noise will now be discussed. For zero mean, temporally 
stationary, and temporally uncprrelated noise the individual frame noise spectrum cross covariance is related to the 
noise power spectral density C NN by 
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{ns(F. tj^(F , 0) = S^^Cnn^F - F ) 

where the individual frame times are displayed explicitly and accommodate the potential of synthesizing images based 
on sets of SpinAp individual frames that may possess common data. 

[0131] The differenced redundancy averaged noise spectrum correlation will now be discussed. For synthesized 
images generated from distinct sets of SpinAp individual frames, the noise spectrum cross-covariance for differenced 
synthesized images is 

. ( N rf (F)N rf (F) N rf {F)N rt (F) \ 

C nn (F)S(F-FK 1- S S fr,,* I £ S 8 k .m 

iN^FjNrflFjj^ i=i .Nrf(FjNrf(F)k-i m=l I 

/N rf {?) N„(f) \ 

= » pfi jEEfr-rf I +2 =-^C™(F)5(F-F) 

N rf (F)N rf (F) I j = i k-r J Nrf(F) 

[01 32J Therefore, the pixel to pixel differenced synthesized noise- for temporally uncorrelated and stationary noise is 

(An syn (R m )An; yn (R)) = 

f d 2 F [ d 2 Fexp[27ii (F R m -F rJ]a(F)A (f )— (F)5(f -F) 
} J NrflF) 

= I d 2 Fexp [2Ki F • (R m - R n )]c nn (F)| A (F)| 2 — 2— 
J NrflFJ 

and, for temporally uncorrelated, temporally stationary, and spatially white, the spatial differenced synthesized noise 
cross-covariance matrix is 

(An syn (R m )AS; yn (R} = J d 2 Fexp [lni F • (R m - R n ))c nn (Fifef -3=r 

; Nrf\FJ . 

= ^_ Sdsef d 2 Fexp (2«i F - (R m - R n )j A(F) | *— 2— . 

4 D x As x D y Ay x J N rf (F) 

[01 33] Embodiments of the present SpinAp moving object detection system consistent with a realization of the meth- 
od illustrated in Figs. 2 and 3 have been reduced to practice and tested. Image data was generated using the embod- 
iment of Fig. 2 to produce a noiseless differenced synthesized moving object signature and a corresponding template 
matched filter output. 

[0134] The image data demonstrated a low signal to noise ratio differenced synthesized image containing a faint 
smear of a moving object, and the output of the -template matched filter, thus demonstrating an excellent output signal 
to noise ratio on the order of 16 to 1 . This demonstrates enhanced receiver operating characteristics that are produced 
using the present invention for this particular class of moving objects. In addition, experimental data has been obtained 
and used to demonstrate the it-wise moving object detection methods described herein. 

[0135] Thus, improved spinning strip (partial) aperture imaging radiometers that synthesize equivalent full circular 
aperture images and detect moving objects from a plurality of rotating strip aperture image measurements, or synthe- 
sized image measurements have been disclosed. It is to be understood that the described 

[01 36] In summary, there is disclosed a spinning aperture imaging radiometer sensor system 1 0 and data processing 
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methods 20 for detecting moving objects derived from a plurality of image frames acquired by a strip aperture imaging 
sensor. A moving object in any invididual image frame results in a motion smear signature in the total synthesized 
image. The motion smear signature is processed to detect the moving objects. One embodiment of the system com- 
prises a rotating strip aperture telescope 11 , a two dimensional detector array 17 that detects images in the telescope's 

5 focal plane, a rotation compensation device 19 that prevents rotational smear during integration time of detectors of 
the array, a signal processor 18 that records a plurality of image frames of a scene imaged by the telescope as it rotates 
around its optical axis, and that implements method(s) 20 for detecting the moving objects present in the recorded 
images. A hierarchy of moving object detection processors 18 and methods 20 is disclosed that includes spatial, tem- 
poral, spatial frequency, and temporal frequency domain detection processors, and is compatible with multi-spectral 

10 background rejection techniques. Selection of the appropriate processing procedure and method depends upon the 
scenario, and the effective signal to noise ratio characteristics of the moving object. The present invention may be 
used for transient event detection. 

[0137] It is to be understood that the described embodiments are merely illustrative of some of the many specific 
embodiments which represent applications of the principles of the present invention. Clearly, numerous and varied 
15 other arrangements may be readily devised by those skilled in the art without departing from the scope of the invention, 
as defined by the appended claims. 

Claims 

20 

1 . A spinning strip aperture radiometer sensor system (10) comprising a telescope (11 ) comprised of a rotating strip 
aperture that rotates around an optical axis and that produces temporally sequential images of an image scene; 

a two-dimensional detector array (17) for detecting images formed at the focal plane of the telescope; and 
25 a signal processor (18) coupled to the detector array for recording a plurality of image frames of the image 

scene as the strip aperture rotates around the optical axis of the telescope, and for synthesizing images from 
the plurality of image frames, 

characterized in that the signal processor (18) is programmed for differencing pairs of images or frames, and 
30 that a matched filter is provided, so as to detect moving objects found in the synthesized images. 

2. The system (10) of claim 1 , characterized in that the signal processor (18) processes at least two synthesized 
images (30) acquired by observing the same scene and detects moving objects found in the synthesized images by: 

35 differencing (40) all possible distinct pairs of synthesized images (30); 

summing (35) the differenced images, and 

processing (50) the resulting differenced and summed data (45) through a matched filter array (55) and a 
comparator (60) whereby a moving object detection is declared when any outputs of the matched filter array 
(55) exceed predetermined thresholds of the comparator. 

40 

3. The system (10) of claim 2, characterized In that the signal processor (18) additionally spatially registers the 
synthesized images due to line of sight errors between measurement sets. 

4. The system (10) of claim 1 , characterized in that the signal processor (18) processes at least two individual line 
45 of sight corrected frames (30-2) resulting from the observation of the same scene and detects moving objects 

found in the images by: 

filtering (31-2) individual frames (30-2) using a system transfer function equalization filter; 
differencing (40-2) distinct pairs of equalized frames; 

50 filtering (41-2) the differenced data using a non-common passband rejection filter; 

summing (35-2) the filtered and differenced data to form a total summed spectrum of filtered differences; 
inverse transforming (42-2) the summed spectra to provide a moving object signature in the spatial domain; and 
processing the resulting summed data (45-2) through a matched filter array (55-2) and comparator (60-2), 
whereby a moving object detection is declared when any outputs of the matched filter array (55-2) exceed 

55 predetermined thresholds of the comparator. 

5. The system (10) of claim 1, characterized in that the signal processor (18) processes a temporal sequence of 
the spatially registered individual image frames (50-3) obtained by observations of the same scene and detects 
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moving objects found in the images by: 

spatially Fourier transforming (51 -3) the spatially registered frames (50-3) to a common spatial frequency grid; 
temporally Fourier transforming (53-3) the spatially registered and Fourier transformed individual frame images 
(52-3) to produce a data set (54-3) associating all temporal frequencies with each spatial frequency; and 
processing (55-3) the temporal and spatial frequency transformed data set (54-3) using a noise optimal 
matched filter detection processor (55-3). 

6. The system (10) of claim 1 , characterized in that the signal processor (18) processes at least two individual line 
of sight corrected frames (30-4) resulting from the observation of the same scene and detects moving objects 
found in the images by: 

filtering (31-4) individual frames (30-4) using a system transfer function equalization filter, 
differencing (40-4) distinct pairs of equalized frames; 

filtering (41-4) the differenced data using a non-common passband rejection filter; 
inverse transforming (42-4) the spectra to provide a moving object signature in the spatial domain; and 
processing (50-4) the resulting individual, equalized, common passband, differenced image data (45-4) 
through a matched filter array (56-4) associated with the set of differenced data (45-4) and then a comparator 
(60-4), whereby a moving object detection is declared when any outputs of the matched filter array (56-4) 
exceed predetermined thresholds of the comparator. 

7. An object detection method (20) for use in an imaging system comprising a spinning strip aperture telescope (11) 
that rotates around an optical axis, a two-dimensional detector array (17) for detecting images located at a focal 
plane of the telescope (11), and a signal processor (18) for processing the detected images, wherein the object 
detection method (20) comprises the following steps: 

recording a plurality of image frames of an image scene imaged by the telescope as the strip aperture rotates 
around the optical axis of the telescope; and 

synthesizing images (30) from the plurality of recorded image frames; 

characterized by differencing pairs of images or frames, and processing the resulting data through a 
matched filter, so as to detect moving objects found in the synthesized images (30). 

8. The method of claim 7, characterized by rotation compensation means (16, 19) for providing a stationary image 
during the integration time of detectors of the detector array (17). 

9. The method of claim 7 or 8, characterized in that the step of detecting moving objects comprises the steps of: 

processing at least two synthesized images (30) acquired by observing the same scene; 
differencing (40) all possible distinct pairs of synthesized images; 
summing (35) the differenced images, and 

processing (50) the resulting differenced and summed data (45) through a matched filter array (55) and then 
a comparator (60), whereby a moving object detection is declared when any outputs of the matched filter array 
(55) exceed predetermined thresholds of the comparator. 

1 0. The method of claim 7 or 8, characterized in that the step of detecting moving objects comprises the steps of: 

processing at least two individual line of sight corrected frames (30-2) resulting from the observation of the 
same scene; 

filtering (31-2) individual frames (30-2) using a system transfer function equalization filter; 
differencing (40-2) distinct pairs of equalized frames; 

filtering (41-2) the differenced data using a non-common passband rejection filter; 
summing (35-2) the filtered and differenced data to form a total summed spectrum of filtered differences; 
inverse transforming (42-2) the summed spectra to provide a moving object signature in the spatial domain; and 
processing (50-2) the resulting summed data (45-2) through a matched filter array (55-2) and comparator 
(60-2), whereby a moving object detection is declared when any outputs of the matched filter array (55-2) 
exceed predetermined thresholds of the comparator. 
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1 1 . The method of claim 7 or 8. characterized in that the step of detecting moving objects comprises the steps of: 

processing at least two individual line of sight corrected frames (30-4) resulting from the observation of the 
s3idg seen©* 

filtering (31-4) individual frames (30-4) using a system transfer function equalization filter; 
differencing (40-4) distinct pairs of equalized frames; 

filtering (41-4) the differenced data using a non-common passband rejection filter; 
inversefransfUng^thesu™^ 

processing (50-1) the resulting individual, equalized, common passband, drfferenced .mage data (£5-4) 
Lough a matched filter array (5fr4) associated with the set of differenced data and then a comparator (60-4). 
whereby a moving object detection is declared when any outputs of the matched filter array (5fr4) exceed 
predetermined thresholds of the comparator. 

12. The method of claim 7 or 8, characterized in that the step of detecting moving objects comprises the steps of: 

processing a temporal sequence of the spatially registered individual image frames (50-3) obtained by the 

observations of the same scene; ... 

spatially Fourier transforming (51 -3) the spatially registered frames (50-3) to a common spafial frequency gnd. 
temporallyFouriertransforming(53-3)thespatially registered and Fourier transformed .nd.yidual frame .mages 
(52-3) to produce a data set (54-3) associating all temporal frequencies with each spat.al frequency; and 
processing (55-3) the temporal and spatial frequency transformed data set (54-3) using a no.se opfimal 
matched filter detection processor. 



Patentanspruche 

1. Radiometersensorsystem (10) mit sich drehender, bandformiger Apertur, mit einem Teleskop 

drehende, bandformige Apertur aufweist, die urn eine optische Achse rofiert und die ze.tl.ch sequenhelle B.lder 
einer Bildszene erzeugt; 

einem zweidimensionalen Detektorarray (17) zum Erfassen von Bildem, die bei der Brennebene desTeleskops 

e^em^nalprozessor (18). der mit dem Detektorarray gekoppelt ist. zum Aufzeichnen einer Vielzahl von 
Bild-Frames der Bildszene, wahrend die bandformige Apertur urn die optische Achse des Teleskops herum 
rotiert. und zum Synthetisieren von Bildem aus der Vielzahl von Bild-Frames. 

dadurch gekennzeichnet, daB der Signalprozessor (18) dazu programmiert ist, Paare von Bildern oder 
Frames zu differenzieren und daft ein angepafttes Filter vorgesehen ist. urn sich bewegende Objekte zu erfassen, 
die sich in den synthetisierten Bildem finden lassen. 

2 System (10) nach Anspruch 1 . dadurch gekennzeichnet. daB der Signalprozessor (18) wenigstens zwei synthe 
tisierte Bilder (30) verarbeitet. die gewonnen werden durch Beobachten der gle.chen Szene. und sich bewegende 
Objekte erfafit, die sich in den synthetisierten Bildem finden lassen, durch: 

Differenzieren (40) samtlicher moglicher verschiedener Paare von synthetisierten Bildem (30); 
Aufsummieren (35) der differenzierten Bilder; und 

Verarbeiten (50) der resultierenden differenzierten und aufsummierten Daten (45) durch e.nen angepaft en 
Filterarray (55) und einen Komparator (60), wobei eine Erfassung eines sich bewegenden Objektes erklart 
wird, wenn irgendwelche Ausgange des angepaftten Filterarrays (55) vorbest.mmte Schwellenwerte des Kom- 
parators Oberschreiten. 

3 System (10) nach Anspruch 2, dadurch gekennzeichnet. daB der Signalprozessor (18) zusatzlich die synthett- 
slerten Bilder raumlich zur Deckung bringt, und zwar auf der Grundlage von Sichtlin.enfehlem zw,schen Satzen 
von Messungen. 

4 System (10) nach Anspruch 1, dadurch gekennzeichnet, daB der Signalprozessor (18) wenigstens zwei indivi- 
" duelle sichtiinien-komgierte Frames (30-2) verarbeitet, die sich aus der Beobachtung der gle.chen Szene ergeben, 

und sich bewegende Objekte, die sich in den Bildem finden lassen, erfafct durch: 
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Fiitem (31-2) von individuellen Frames (30-2) und unter Verwendung eines Systemubertragungsfunktion-Aus- 
gleichsfilters; 

Differenzieren (40-2) von verschiedenen Paaren von ausgeglichenen Frames; 

Fiitem (41-2) der differenzierten Daten unter Verwendung eines Zuruckweisungsftlters fur ein nicht gemein- 
sames DurchlaRband; 

Aufsummieren (35-2) der gefilterten und differenzierten Daten, urn ein aufsummiertes Gesamtspektrum von 
gefilterten Differenzen zu bilden; 

Umkehrtransformieren (42-2) der aufsummierten Spektren, urn eine Signatur eines sich bewegenden Objektes 
in der Raumebene bereitzustellen; und 

Verarbeiten der resultierenden, aufsummierten Daten (45-2) mitteis eines angepafiten Filterarrays (55-2) und 
eines Komparators (60-2), wobei eine Erfassung eines sich bewegenden Objektes erklart wird, wenn irgend- 
welche Ausgange des angepafiten Filterarrays (55-2) vorbestimmte Schwellenwerte des Komparators uber- 
schreiten. 

System (10) nach Anspruch 1, dadurch gekennzeichnet, daB der Signalprozessor (18) eine zeitliche Folge der 
raumlich zur Deckung gebrachten, individuellen Bild-Frames (50-3) verarbeitet, die erhalten werden durch Beob- 
achtungen der gleichen Szene, und sich bewegende Objekte, die sich in den Bildem finden lassen, erfafit durch: 

raumliches Fourier-Transformieren (51-3) der raumlich zur Deckung gebrachten Frames (50-3), auf ein ge- 
meinsames raumliches Frequenzgitter; 

zeitliches Fourier-Transformieren (53-3) der raumlich zur Deckung gebrachten und fourier-transformierten 
individuellen Frame-Bilder (52-3), urn einen Datensatz (54-3) zu erzeugen, der samtliche zeitlichen Frequen- 
zen jeder raumlichen Frequenz zuordnet; und 

Verarbeiten (55-3) des zeitlich und raumlich frequenztransformierten Datensatzes (54-3) unter Verwendung 
eines rausch-optimalen angepafiten Filtererfassungsprozessors (55-3). 

System (10) nach Anspruch 1, dadurch gekennzeichnet, dad der Signalprozessor (18) wenigstens zwei indivi- 
duelle sichtlinien-korrigierte Frames (30-4) verarbeitet, die sich ergeben aus der Beobachtung der gleichen Szene, 
und sich bewegende Objekte erfafit, die sich in den Bildem finden lassen, durch: 

Fiitem (31-4) von individuellen Frames (30-4) unter Verwendung eines Systemubertragungsfunktion-Aus- 
gleichsfilters; 

Differenzieren (40-4) von verschiedenen Paaren von ausgeglichenen Frames; 

Fiitem (41-4) der differenzierten Daten unter Verwendung. eines Zurtickweisungsfi Iters fur ein nicht gemein- 
sames Durchlafiband; 

Umkehrtransformieren (42-4) der Spektren, urn eine Signatur eines sich bewegenden Objektes in der Raum- 
ebene bereitzustellen; und 

Verarbeiten (50-4) der sich ergebenden individuellen, ausgeglichenen, ein gemeinsames Durchlafiband auf- 
weisenden, differenzierten Bilddaten (45-4) mitteis eines angepaBten Filterarrays (56-4), der dem Satz von 
differenzierten Daten (45-4) zugeordnet ist, und dann mitteis eines Komparators (60-4), wobei eine Erfassung 
eines sich bewegenden Objektes erklart wird, wenn irgendwelche Ausgange des angepafiten Filterarrays 
(56-4) vorbestimmte Schwellenwerte des Komparators uberschreiten. 

Verfahren (20) zur Erfassung eines Objektes, zur Verwendung in einem bildgebenden System, das ein Teleskop 
(11 ) mit einer sich drehenden, bandformigen Apertur aufweist, das urn eine optische Achse rotiert, das femer einen 
zweidimensionalen Detektorarray (17) aufweist, zum Erfassen von Bildem, die bei der Brennebene des Teleskopes 
(11) vorliegen, und das ferner einen Signalprozessor (18) aufweist, zum Verarbeiten der erfafcten Bilder, wobei 
das Verfahren (20) zur Erfassung eines Objektes die folgenden Schritte aufweist: 

Aufzeichnen einer Vielzahl von Bild-Frames einer Bildszene, die durch das Teleskop abgebildet werden, wenn 
die bandformige Apertur um die optische Achse des Teleskops rotiert; und 
Synthetisieren von Bildem (30) aus der Vielzahl von aufgezeichneten Bild-Frames; 

gekennzeichnet durch Differenzieren von Paaren von Bildern oder Frames und durch Verarbeiten der sich 
ergebenden Daten mitteis eines angepaBten Filters, um sich bewegende Objekte zu erfassen, die sich in den 
synthetisierten Bildem (30) finden lassen. 

Verfahren nach Anspruch 7, gekennzeichnet durch Drehkompensationsmittel (16, 19) zum Bereitstellen eines 
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stationaren Bildes wahrend der Integrationszeit von Detektoren des Detektorarrays (17). 

9. Verfahren nach Anspruch 7 oder 8, dadurch gekennzeichnet, daB der Schritt des Erfassens von sich bewegenden 
Objekten die Schritte aufweist: 

5 

Verarbeiten von wenigsten zwei synthetisierten Bildem (30), die gewonnen werden durch Beobachten der 
gleichen Szene; 

Differenzieren (40) samtlicher moglicher verschiedener Paare von synthetisierten Bildem; 
Aufsummieren (35) der differenzierten Bilder; und 
w Verarbeiten (50) der sich ergebenden, differenzierten und aufsummierten Daten (45) mittels eines angepaBten 

Filterarrays (55) und dann mittels eines Komparators (60), wobei eine Erfassung eines sich bewegenden 
Objektes erklart wird, wenn irgendwelche Ausgange des angepaBten Filterarrays (55) vorbestimmte Schwel- 
lenwerte des Komparators uberschreiten. 

15 10. Verfahren nach Anspruch 7 Oder 8, dadurch gekennzeichnet, daB der Schritt des Erfassens von sich bewegenden 
Objekten die Schritte aufweist: 

Verarbeiten von wenigstens zwei individuellen sichtlinienkorrigierten Frames (30-2), die sich aus der Beob- 
achtung der gleichen Szene ergeben; 
20 Filtem (31-2) von individuellen Frames (30-2) unter Verwendung eines Systemubertragungsfunktion-Aus- 

gleichsfilters; 

Differenzieren (40-2) von verschiedenen Paaren von ausgeglichenen Frames; 

Filtem (41-2) der differenzierten Daten unter Verwendung eines Zuruckweisungsfilters fur ein nicht gemein- 
sames DurchlaBband; 

25 Aufsummieren (35-2) der gefilterten und differenzierten Daten, um ein aufsummiertes Gesamtspektrum von 

gefilterten Differenzen zu bilden; 

Umkehrtransformieren (42-2) der aufsummierten Spektren, um eine Signatur eines sich bewegenden Objektes 
in der Raumebene bereitzustellen; und 

Verarbeiten (50-2) der sich ergebenden aufsummierten Daten (45-2) mittels eines angepaBten Filterarrays 
30 (55-2) und eines Komparators (60-2), wobei eine Erfassung eines sich bewegenden Objektes erklart wird, 

wenn irgendwelche Ausgange des angepaBten Filterarrays (55-2) vorbestimmte Schwellenwerte des Kom- 
parators uberschreiten. 

1 1 . Verfahren nach Anspruch 7 oder 8, dadurch gekennzeichnet daB der Schritt des Erfassens von sich bewegenden 
35 Objekten die Schritte aufweist: 

Verarbeiten von wenigstens zwei individuellen, sichtltnien-korrigierten Frames (30-4), die sich ergeben aus 
der Beobachtung der gleichen Szene; 

Filtem (31-4) von individuellen Frames (30-4) unter Verwendung eines System ubertragu ngsfunktion-Aus- 
40 gleichsfi Iters; 

Differenzieren (40-4) von verschiedenen Paaren von ausgeglichenen Frames; 

Filtem (41-4) der differenzierten Daten unter Verwendung eines Zuruckweisungsfilters fur ein nicht gemein- 
sames DurchlaBband; 

Umkehrtransformieren (42-4) der aufsummierten Spektren, um eine Signatur eines sich bewegenden Objektes 

45 in der Raumebene bereitzustellen; und 

Verarbeiten (50-4) der sich ergebenden individuellen, ausgeglichenen, ein gemeinsames DurchlaBband auf- 
weisenden differenzierten Bilddaten (45-4) mittels eines angepaBten Filterarrays (56-4), der dem Satz von 
differenzierten Daten zugeordnet ist, und dann mittels eines Komparators (60-4), wobei eine Erfassung eines 
sich bewegenden Objektes erklart wird, wenn irgendwelche Ausgange des angepaBten Filterarrays (56-4) 

50 vorbestimmte Schwellenwerte des Komparators uberschreiten. 

12. Verfahren nach Anspruch 7 Oder 8, dadurch gekennzeichnet, daB der Schritt des Erfassens von sich bewegenden 
Objekten die Schritte aufweist: 

55 Verarbeiten einer zeitlichen Folge der raumlich zur Dekkung gebrachten, individuellen Bild-Frames (50-3), die 

erhalten werden durch die Beobachtungen der gleichen Szene; 

raumliches Fourier-Transformieren (51-3) der raumlich zur Deckung gebrachten Frames (50-3), auf ein ge- 
meinsames raumliches Frequenzgitter; 
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zeitliches Fourier-Transformieren (53-3) der raumlich zur Deckung gebrachten und fourier-transfomierten in- 
dividuellen Frame-Bilder (52-3), urn einen Datensatz (54-3) zu erzeugen, der samtliche zeitlichen Frequenzen 
jeder raumlichen Frequenz zuordnet; und 

Verarbeiten (55-3) des zeitlich und raumlich frequenztransformierten Datensatzes (54-3) unter Verwendung 
5 eines rausch-optimalen angepaRten Filtererfassungsprozessors. 



Revendications 

10 1 . Systeme (1 0) de capteur de radiometre a ouverture tournante en forme de bande comprenant un telescope (11) 
constitue d'une ouverture tournante en forme de bande qui toume autour d'un axe optique et qui produit tempo- 
rellement des images sequentielles d'une scene d'image ; 

un groupement bidimensionnel (17) de detecteurs destine a detecter des images formees au niveau du plan 
15 focal du telescope ; et 

un processeur (18) de signal associe au groupement de detecteurs pour enregistrer une pluralite de vues 
d'image de la scene d'image lorsque I'ouverture en forme de bande toume autour de I'axe optique du telescope, 
et a synthetiser des images a partir de la pluralite de vues d'image ; 

20 caracterise en ce que le processeur (18) de signal est programme pour differencier des paires d'images 

ou de vues, et en ce qu'il est prevu un filtre adapte, de facon a detecter des objets en mouvement trouves dans 
les images synthetisees. 

2. Systeme (10) selon la revendication 1, caracterise en ce que le processeur (18) de signal traite au moins deux 
25 images synthetisees (30) acquises par I'observation de la meme scene et detecte des objets en mouvement trouves 

dans les images synthetisees : 

en differencial (40) toutes les paires distinctes possibles d'images synthetisees (30) ; 
en totalisant (35) les images differenciees ; et 
30 en traitant (50) les donnees resultantes (45) differenciees et totalisees au moyen d'un groupement (55) de 

filtres adaptes et d'un comparateur (60), en declarant ainsi une detection d'objet en mouvement lorsque des 
sorties quelconques du groupement (55) de filtres adaptes depassent des seuils predetermines du compara- 
teur. 

35 3. Systeme (10) selon la revendication 2, caracterise en ce que le processeur (18) de signal enregistre spatialement, 
en outre, les images synthetisees dues a des erreurs de ligne de visee entre des ensembles de mesures. 

4. Systeme (10) selon la revendication 1, caracterise en ce que le processeur (18) de signal traite au moins deux 
vues (30-2) individuelles corrigees en ce qui concerne la ligne de visee, resultant de Pobservation de la meme 
scene, et detecte des objets en mouvement trouves dans les images : 

en filtrant (31-2) les vues individuelles (30-2) en utilisant un filtre d'egalisation de fonction de transfert du 
systeme ; 

en differencial (40-2) des paires distinctes de vues egalisees ; 
45 en filtrant (41-2) les donnees differenciees en utilisant un filtre a rejection de bande passante non commune ; 

en totalisant (35-2) les donnees filtrees et differenciees pour former un spectre totalise global des differences 
filtrees ; 

en transformant de maniere inverse (42-2) les spectres totalises pour foumir une signature d'objet en mou- 
vement dans le domaine spatial ; et 
50 en traitant les donnees totalisees resultantes (45-2) au moyen d'un groupement (55-2) de filtres adaptes et 

d'un comparateur (60-2), en declarant ainsi une detection d'objet en mouvement lorsque des sorties quelcon- 
ques du groupement (55-2) de filtres adaptes depassent des seuils predetermines du comparateur. 

5. Systeme (10) selon la revendication. 1 , caracterise en ce que le processeur (18) de signal traite une sequence 
55 temporelle des vues (50-3) d'image individuelles enregistrees spatialement obtenues par des observations de la 

meme scene et detecte des objets en mouvement trouves dans les images : 

en transformant spatialement (51-3) par transformee de Fourier les vues enregistrees spatialement (50-3) 
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dans une grille de frequences spatiales commune ; 

en transformant tempore! lement (53-3) par transformee de Fourier les Images (52-3) de vues individuelles 
enreglstres spatialement et transformees par transformee de Fourier pour produire un ensemble (54-3) de 
donnees associant toutes les frequences temporelles avec chaque frequence spatiale ; et 
5 en traitant (55-3) I'ensemble (54-3) de donnees transformees a frequences temporelles et spatiales en utilisant 

un processeur (55-3) de detection de bruit a filtre adapte optimal. 

6. Systeme (10) selon la revendication 1, caracterise en ce que le processeur (18) de signal traite au moins deux 
vues (30-4) individuelles corrigees en ce qui conceme la ligne de visee, resultant de I'observation de la meme 

10 scene, et detecte des objets en mouvement trouves dans les images : 

en filtrant (31-4) les vues individuelles (30-4) en utilisant un filtre d'egalisation de fonction de transfer! du 
systeme ; 

en differencial (40-4) des paires distinctes de vues egalisees ; 
15 en filtrant (41-4) les donnees differenciees en utilisant un filtre a rejection de bande passante non commune ; 

en transformant de maniere inverse (42-4) les spectres pour foumir une signature d'objet en mouvement dans 
le domaine spatial ; et 

en traitant (50-4) les donnees (45-4) individuelles resultantes, egalisees, differenciees, a bande passante 
commune, au moyen d'un groupement (56-4) de filtres adaptes associe a I'ensemble de donnees (45-4) dif- 
20 ferenciees, et au moyen, ensuite, d'un comparateur (60-4), en declarant ainsi une detection d'objet en mou- 

vement lorsque des sorties quelconques du groupement (56-4) de filtres adaptes depassent des seuils pre- 
determines du comparateur. 

7. Precede (20) de detection d'objets pour utilisation dans un systeme de saisie d'image comprenant un telescope 
25 (1 1 ) a ouverture toumante en forme de bande qui toume autour d'un axe optique, un groupement bidimensionnel 

(17) de detecteurs destine a detecter des images situees au niveau d'un plan focal du telescope (11), et un pro- 
cesseur (18) de signal destine a traiter les images detectees, dans lequel le precede (20) de detection d'objets 
comprend les etapes consistant : 

30 a enregistrer une pluralite de vues d'image d'une scene d'image dont I'image est saisie par le telescope lorsque 

I'ouverture en forme de bande toume autour de I'axe optique du telescope ; et 
a synthetiser des images (30) a partir de la pluralite de vues d'image enregistree ; 

caracterise par le fait de differencier des paires d'images ou de vues, et de traiter les donnees resultantes 
35 au moyen d'un filtre adapte, de facon a detecter des objets en mouvement trouves dans les images synthetisees 

(30). 

8. Precede selon la revendication 7, caracterise par un moyen (1 6, 1 9) de compensation de rotation destine a foumir 
une image fixe pendant le temps d'integration des detecteurs du groupement (17) de detecteurs. 

40 

9. Precede selon la revendication 7 ou 8, 

caracterise en ce que I'etape de detection d'objets en mouvement comprend les etapes consistant : 

a traiter au moins deux images synthetisees (30) acquises par I'observation de la meme scene ; 
45 a differencier (40) toutes les paires distinctes possibles d'images synthetisees ; 

a total iser (35) les images differenciees ; et 

a traiter (50) les donnees resultantes (45) differenciees et totalisees au moyen d'un groupement (55) de filtres 
adaptes et au moyen, ensuite, d'un comparateur (60), en declarant ainsi une detection d'objet en mouvement 
lorsque des sorties quelconques du groupement (55) de filtres adaptes depassent des seuils predetermines 
50 du comparateur. 

10. Precede selon la revendication 7 ou 8, caracterise en ce que I'etape de detection d'objets en mouvement com- 
prend les etapes consistant : 

a traiter au moins deux vues (30-2) individuelles corrigees en ce qui conceme la ligne de visee, resultant de 
55 Pobservation de la meme scene : 

a filtrer (31-2) les vues individuelles (30-2) en utilisant un filtre legalisation de fonction de transfert du systeme ; 
a differencier (40-2) des paires distinctes de vues egalisees ; 



38 



EP 0 802 426 B1 

a filtrer (41-2) les donnees differences en utilisant un filtre a rejection de bande passante non commune ; 
a totaliser (35-2) les donnees filtrees et differenciees pour former un spectre totalise global des differences 
filtrees ; 

a transformer de maniere inverse (42-2) les spectres totalises pour foumir une signature d'objet en mouvement 
5 dans le domaine spatial ; et 

a traiter (50-2) les donnees total isees resultantes (45-2) au moyen d'un groupement (55-2) de filtres adaptes 
et d'un comparateur (60-2), en declarant ainsi une detection d'objet en mouvement lorsque des sorties quel- 
conques du groupement (55-2) de filtres adaptes depassent des seuits predetermines du comparateur. 

10 11. Procede selon ta revendication 7 ou 8, caracterise en ce que Petape de detection d'objets en mouvement com- 
prend les etapes consistant : 

a traiter au moins deux vues (30-4) individuelles corrigees en ce qui conceme la ligne de visee, resultant de 
('observation de la meme scene : 

15 a filtrer (31-4) les vues individuelles (30-4) en utilisant un filtre d'egalisation de fonction detransfertdu systeme ; 

a differencier (40-4) des paires distinctes de vues egalisees ; 

a filtrer (41-4) les donnees differenciees en utilisant un filtre a rejection de bande passante non commune ; 
a transformer de facon inverse (42-4) les spectres totalises pour donner une signature d'objet en mouvement 
dans le domaine spatial ; et 

20 a traiter (50-4) les donnees d'image (45-4) differenciees individuelles resultantes, egalisees, a bande passante 

commune, au moyen d'un groupement (56-4) de filtres adaptes associe a I'ensemble des donnees differen- 
ciees et au moyen, ensuite, d'un comparateur (60-4), en declarant ainsi une detection d'objet en mouvement 
lorsque des sorties quelconques du groupement (56-4) de filtres adaptes depassent des seuils predetermines 
du comparateur. 

25 

12. Procede selon la revendication 7 ou 8, caracterise en ce que I'etape de detection d'objets en mouvement com- 
prend les etapes consistant : 

a traiter une sequence temporetle des vues (50-3) d'image individuelles enregistrees spatialement, obtenues 
30 par les observations de la meme scene ; 

a transformer spatialement (51 -3) par transformee de Fourier, les vues enregistrees spatialement (50-3) en 
une grille de frequences spatiales commune ; 

a transformer temporellement (53-3) par transformee de Fourier les images (52-3) de vue individuelles enre- 
gistrees spatialement et transformees par transformee de Fourier, pour produire un ensemble (54-3) de don- 
35 nees associant toutes les frequences temporelles a chaque frequence spatiale ; et 

a traiter (55-3) I'ensemble (54-3) de donnees transformees de frequences temporelles et spatiales, en utilisant 
un processeur de detection par filtre adapte a bruit optimal. 
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